Asymptotics of the ground state energy of 
(N heavy molecules and related topics 

T— I 

o 

^ Victor Ivrii 

O October 5, 2012 

4^ Abstract 
Oh 

We consider asymptotics of the ground state energy of heavy 
atoms and molecules in the strong external magnetic field and derive 
^ it including Schwinger and Dirac corrections (if magnetic field is 

not too strong). In the next instalment we extend this paper to 
^ consider also related topics: an excessive negative charge, ionization 

energy and excessive negative charge when atoms can still bind into 
molecules. 
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Introduction 

In this paper we repeat analysis of the previous paper [IvrlO]^^ but in the 
case of the constant external magnetic field^-*. 

. 1 Framework 

Let us consider the following operator (quantum Hamiltonian) 
(0.1.1) H = H^ := J2 Ha.v..,+ Y1 l^i-^'^l"' 

l<j<N l<J<k<N 

on 

(0.1.2) Sj= /\ Ji^, = ^^R^'X") 

l<n<N 

with 

(0.1.3) Hv.A = ((/V - A) ■ cr)' - V{x) 

describing N same type particles in the external field with the scalar potential 
— V and vector potential A[x), and repulsing one another according to the 
Coulomb law. 

Here xj e R'' and (xi,...,xa/) G M""^, potentials V{x) and A{x) are 
assumed to be real-valued. Except when specifically mentioned we assume 
that 

(0.1.4) V{x)= J2 

l<k<M 



\X 



where > and yk are charges and locations of nuclei. Here cr = 
(ai, 0'2, cTd), cT/c are q x g-Pauli matrices. 

So far in comparison with the previous paper [IvrlO] we only changed 
(24.1.3) of [Ivrll] to (0.1.3) introducing magnetic field. Now spin enters not 
only in the definition of the space but also into operator through matrices 



^) Coinciding with Chapter 24 of of [Ivrll]. 

^•^ ActuaUy we need a magnetic field either sufficiently weak or close to a constant on 
the very small scale. 
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(Jk. Since we need of = 3 Pauli matrices it is sufficient to consider q = 2 but 
we will consider more general case as well (but q should be even). 

Remark 0.1.1. In the case of the the constant magnetic field V x A 
(0.1.5) Ha,v = (-/V - A{x)f + a-W xA - V{x) 

In the case d = 2 this operator downgrades to 

(0.1.6) H^y = (-'V - A{x)f + (73(V xA)- V{x) 

Again, let us assume that 
(0.1.7) Operator H is self-adjoint on S). 

As usual we will never discuss this assumption. 

0.2 Problems to consider 

As in the previous Chapter we are interested in the ground state energy 
E = Ea/ of our system i.e. in the lowest eigenvalue of the operator H = H/v 
on S^: 

(0.2.1) E := infSpecH on Sj; 

more precisely we are interested in the asymptotics of E/v = E(y; Z; N) as V 
is defined by (0.1.4) and A/ x Z := Zi + Z2 + ... + Zm 00 and we are going 
to prove that^^ E is equal to Magnetic Thomas-Fermi energy £q^, possibly 
with Scott and Dirac-Schwinger corrections and with appropriate error. 
We are also interested in the asymptotics for the ionization energy 

(0.2.2) l/v := -E/v + Ea/-i 

and we also would like to estimate maximal excessive negative charge 
(0.2.3) max N - Z. 

N: l„>0 

Under reasonable assumption to the minimal distance between nuclei. 
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All these questions so far were considered in the framework of the fixed 
positions yi, ... , yw but we can also consider 

(0.2.4) E = E/v = E(y; Z;N) = E+ U{r, Z) 

with 

(0.2.5) 
and 

(0.2.6) E(Z;/V)= inf E{y; Z; N) 

yi vm - 

and replace l/v by //y = — E/y + E/y-i and modify all our questions accord- 
ingly. We call these frameworks fixed nuclei model and free nuclei model 
respectively. 

In the free nuclei model we can consider two other problems: 

- Estimate from below minimal distance between nuclei i.e. 

min |y^-y„/| 

l<m<m'<M 

for which such minimum is achieved; 

- Estimate maximal excessive positive charge 

(0.2.7) max{Z - /V : E < min ^ E(Z„; A/„)} 

A/i,...,Wm: l<m<M 
Ni+...Nm=N 

for which molecule does not disintegrates into atoms. 



0.3 Magnetic Thomas-Fermi theory 

As in the previous Chapter the first approximation is the Hartrcc-Fock (or 
Thomas-Fcrmi) theory. Let us introduce the spacial density of the particle 
with the state e S^: 

(0.3.1) p{x) = p^{x) = N J |M/(x, X2, ... , x/v)|' dx2-- dx^. 
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Let us write the Hamiltonian, describing the corresponding "quantum hq- 
uid": 

(0.3.2) Sb{p) = j Ts{p{x)) dx-j V{x)p{x) dx + ^D(p, p), 

with 

(0.3.3) D(p, p) = JJ\x- y\-'p{x)p{y) dxdy 

where Tg is the energy density of a gas of noninter acting electrons: 
(0.3.4) tb{p) = sup{pw - P{w)) 

w>0 

is the Legendre transform of the pressure Pb{w) given by the formula 
(0.3.5) Pb{w) = + ^(1/1/ - 2ie)f) 

with xi = {27r)~^q, {37r^)~^q for d = 2,3 respectively. 

The classical sense of the second and the third terms in the right-hand 
expression of (0.3.2) is clear and the density of the kinetic energy is given by 
tb{p) in the semiclassical approximation (see remark 0.3.1). So, the problem 
is 

(0.3.6) Minimize functional Sb{p) defined by (0.3.2) under restrictions: 

(0.3.7)^2 P^O' Jpdx<N. 

The solution if exists is unique because functional £b{p) is strictly convex 
(see below). The existence and the property of this solution denoted further 
by pg"^ is known in the series of physically important cases. 

Remark 0.3.1. \i w is the negative potential then 
(0.3.8) tre(x,x,0) w Pe(>^) 

defines the density of all non-interacting particles with negative energies at 
point X and 

(0.3.9) I Tdrtre{x,x,T)dx - f PB{w)dx 

J — oo J 

is the total energy of these particles; here ~ means "in the semiclassical 
approximation" . 
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We consider in the case of c/ = 3 a large (heavy) molecule with potential 
(24.1.4) of [Ivrll]. It is well-known^) that 

Proposition 0.3.2. (i) For V{x) given by (0.1.4) minimization problem 
(0.3.6) has a unique solution p = p^^ ; then denote E^q := Eb{,Pb^', 

(ii) Equality in (0.3.7)2 holds if and only if N < Z .= Z^; 

(Hi) Further, p^^ does not depend on N as N > Z; 

(iv) Thus 



0.4 Main results sketched and plan of the 
chapter 

In the first half of this Chapter we derive asymptotics for ground state energy 
and justify Thomas- Fermi theory. As construction of Section 24.2 of [Ivrll] 
works with minimal modifications (see Section 5) in the magnetic case as 
well we start immediately from magnetic Thomas- Fermi theory in Section 1. 
We discover that there are three different weak magnetic field case 

6 <^ Zs when x Zs and S^^ = ^^0^*^(1 + o(l)), a strong magnetic field 
case e > Z5 when S^^ x bIzI and S^^ = £^{1 + o(l)) where is 
Thomas-Fermi potential derived as Pb{w) = (cf. (0.3.5)), and an 

4 

intermediate case B ^ Zi. 

Then we apply semiclassical methods (like in Section 24.4) of [Ivrll] 
albeit now analysis is way more complicated due to two factors: semi- 
classical theory of magnetic Schrodinger operator is more difficult than 
the corresponding theory for non-magnetic Schrodinger operator and also 
Thomas-Fermi potential W'^^ is not very smooth in the magnetic case, so 
we need to approximate it by a smooth one (on a microscale). We discover 
that both semiclassical methods and Thomas-fermi theory are relevant only 
as B <^ Z^. Case of the superstrong magnetic field B ^ Z^ was considered 
in E. H. Lieb, J. P. Solovej and J. Yngvarsson [LSYl] and we hope to cover 
it in the Chapter 27"'" 

4) Section IV of E. H. Lieb, J. P. Solovej and J. Yngvarsson [LSY2]. 



(0.3.10) 




l<m<M 
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^ Pending 



First of all, in Section 2 we consider the case M = 1; then Thomas- Fermi 
potential Wj^ is non-degenerate and in this case we derive sharp spectral 
asymptotics. 

Next, in Section 3 we consider the case M > 2 but only in the zone 
{W^^ + iy>B} where u is a chemical potential and B is an intensity of the 
magnetic field. A certain weaker non- degeneracy condition is satisfied due 
to Thomas-Fermi equation and we derive almost sharp spectral asymptotics. 

Then in Section 4 we consider the case M > 2 in the zone { W'^^ + < 6} 
containing the boundary of supp pg"^; this is the most difficult case to analyze 

2 

and our remainder estimates are not sharp unless N > Z — CZi. 

Finally, in Section 5 we derive asymptotics of the ground state energy. 
Their precision (or lack of it) follows from the precision of the corresponding 
semiclassical results; so our results in the case M = 1 are sharp, but results 

2 

in the case M >2 (especially if A/ < Z — CZa) are not. 

In the second half^ of this Chapter we consider related problems. In 
Section 6 (cf. Section 24.5) of [Ivrll] we consider negatively charged systems 
(A/ > Z) and estimate both ionization energy l/v and excessive negative 
charge {N — Z)+^\ 

In Section 7 (cf. Section 24.6) of [Ivrll] we consider positively charged 
systems (A/ < Z) and estimate the remainder |Ia/ + z^| in the formula l/v ~ —i'] 
as M > 2 we also estimate excessive positive charge (Z — A/)+ when atoms 
can be bound into molecule^^. 

We will assume that 



(0.4.1) 



,-1 



N <Z^<cN 



Vm 



M. 



1 Magnetic Thomas- Fermi theory 
1.1 Framework and existence 

The Thomas- Fermi theory is well developed in the magnetic case as well albeit 
in the lesser degree than in the non-magnetic one. The most important source 
now is Section IV of E. H. Lieb, J. P. Solovej and J. Yngvarsson [LSY2]. 

Again as in the previous Paper [IvrlO] to get the best lower estimate 
for the ground state energy (neglecting semiclassical errors) one needs to 

In (magnetic) Thomas-Fermi theory both answers are 0. 
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maximize functional ^b,*{^ + ^) defined by (24.3.1) of [Ivrll] albeit with 
the pressure Pb{w) given for d = 2,3 by (0.3.5). Formulae (24.3.2) and 
(24.3.3) of [Ivrll] also remain valid. 

Further, to get the best upper estimate (neglecting semiclassical errors) 
one needs to minimize functional ^*b{p', i^) defined by (24.3.4) of [Ivrll] 
where (24.3.4) of [Ivrll] remains valid with P replaced by Pg and respectively 
r(p') replaced by tb{p') which is Legendre transformation of Pg (see (0.3.4)). 

Since Pg is given by much more complicated expression (0.3.5) rather 
than (24.3.6) of [Ivrll], and respectively 

(1.1.1) P'siw) = ^^,B(^wf' + ^(w - 2jB)i-') 

(cf. (24.3.6)2) [Ivrll] there is no explicit expression for rg similar to 
(24.3.7) of [Ivrll]. 

Remark 1.1.1. (i) e(x) = |V x A[x)\] 

(ii) From now on we will assume that d = 3; 

(iii) Pg is a strictly convex function and therefore rg is also a strictly convex 
function^-*; 

(iv) Pg(iv) Pq{w), P'b{w) Po{w) and rg(p) ro(p) as 6 where 
(without subscript "0") the limit functions have been defined by (24.3.6) 
and (24.3.7) of [Ivrll] respectively. 

Remark 1.1.2. (i) Alternatively we minimize £b{p) = 'i'g(P'O) under as- 
sumptions 

(1.1.2) ^ 2 P > 0' J pdx< N; 

(ii) So far in comparison with the previous paper [IvrlO] we changed only 
definition of Pg(i/i/) and rg(p) respectively. Note that Pg(i/i/) belongs to ^2 
(as d = 2,3) as function of w; this statement will be quantified later; 

As d — 2, Pb is a convex piecewise linear function and therefore rg is also a convex 
function. 
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(iii) While not affecting existence (witli equality in (1.1.2)^^ iff N < Z) and 

4 

uniqueness of solution, it affects other properties, especially as 6 > Zs. 

Proposition 1.1.3. In our assumptions for any fixed u < statements 
(i)-(vii) of proposition 24-3.1 of [Ivrll] hold. 

Proof. The proof is the same as of proposition 24.3.1 of [Ivrll]. The proof 
that threshold u = matches to A/ = Z are theorems 4.9 and 4.10 of Section 
IV of E. H. Lieb, J. P. Solovej and J. Yngvarsson [LSY2]. □ 

Note that (24.3.8)-(24.3.9) and (24.3.10) of [Ivrll] become 

(1.1.3) p= —A{W - V) = P'b{W + 

An 

(1.1.4) 1/1/ = o(l) as |x| oo 
and 

(1.1.5) M{v) = j P'b{W + v) dx 
respectively. 

Similarly, proposition 24.3.2 of [Ivrll] remains true: 

Proposition 1.1.4. For arbitrary W the following estimates hold with ab- 
solute constants cq > and Cq: 

(1.1.6) toDi^p - p"^ , p - p"^) < cDg,,(l/l/TF + ^,)_(Dg ^(i/i/ + ^,) < 

QD{p-p',p-p') 

and 

(1.1.7) eoD(p' - p' - p-"') < %{p, u) - ct)*e(pT^ u) < 

CoD{p-p',p-p') 

with p=lA{W- V), p' = P'b{W + u). 

Proof. This proof is rather obvious as well. □ 
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1.2 Properties 

Proposition 1.2.1. The solution of the magnetic Thomas-Fermi problem 
has the following scaling properties 

(1.2.1) W'^^[x; Z; y; B; N; q) = 

qlNlw^^{q-^N'^x; N'^Z; qi/V^y; q-^N-^B; 1;1), 

(1.2.2) p'^^{x; Z; y; B; N; q) = 

q^N'^p^^q-^N-^x; N-^Z; qi/V^y; q-^N-lB; 1;1), 

(1.2.3) £^^{Z; y; B; N; q) = ql Nh'^^ {N'^ Z; qi/V^y; q-^N-lB; 1; 1), 

(1.2.4) z/^^(Z; y; 6; A/; g) = q'^ N'^ u'^^ {N'^ Z; c/i/V^y; q-'^N-'^B; 1; 1) 

where v^^ = v is the chemical potential; recall that Z = (Zi, ... , Z/v?) arac? 
y = (yi, ... , y/v?) are arrays and parameter q also enters into Thomas-Fermi 
theory. 

In particular, v^^ and B scale the same way. 
Proof. Proof is trivial by scaling. □ 

4 4 

Now one can guess that there are two cases 6 <^ Za and B ^ Z^ (recall 
that A/ X Z) in which magnetic Thomas-Fermi theory looks very different 

4 _ 

(and also an intermediate case 6 ~ Zs). To explain this difference let us 
consider one atom case: 

First of all recall that if 6 = and 
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respectively where the first factors are 
spacial densities of the charge and (negative) Thomas-Fermi energy respec- 
tively and therefore zone ^ x Z~3 provides the main contributions into 
both. 

Therefore, if in this main zone 6 <^ 1/1/"" x Zs we guess that the 
magnetic theory is similar to non-magnetic one. and actually it is true. 

However, let us study an atomic case rigorously. Let M = 1, = and 
N <Z. Then 

(1.2.5) is a spherically symmetric, and it is monotone non-increasing 
function of |x|; Wj^ — >■ -|-0 as |x| — >■ oo; 

(1.2.6) Wl^{x) <-v ^ Wl^ = \x\-\Z - N). 

Really, (1.2.5) is obvious and (1.2.6) follows from it and Newton screening 
theorem. 

4 4 

Two propositions below treat cases 6 < Zs and 6 > Zs respectively; in 

4 4 

the former case there is another fork: 6 < (Z — N)^ and 6 > (Z — N)^. 
Proposition 1.2.2. Let M = 1, y„ = 0, A/ x Z„ and S < Zt . 

(i) Then 

(1.2.7) Wj^ < min(Z|x|-\ C|x|-^) 
and 

(1.2.8) pY < Cmin(zi|x|-i + eZ5|x|-i, \x\~^ + B\x\-^); 
(a) There exists 

(1.2.9) r„xmin(6-5,(Z-/V)7) 

such that ^ —u as |x| ^ and then p^J = iff x > r^; 

(Hi) (1.2.7) and (1.2.8) become equivalencies (xj as \x\ < (1 — e)rm; 

4 _ 1 4 

(iv) B <{Z - N)l implies rrr, >^ {Z - /\/)+% z/ x (Z - N)l and 

(1.2.10) l/l/"^'= + i/x(Z-/V)|.(r„-|x|), 

- a,,, W^^ X (Z - N)l as (1 - e)r„ < |x| < 
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(v) B>{Z- N)l implies x 6"^, i/ x (Z - N)+B'i < B and 

(1.2.11) 1/1/"^'= + ^/ X e2(r„ - |x|)^ + ei(Z - /V)+(r„ - |x|) 

- 1/1/"^'^ X B^ir^r, - \x\f + fii(Z - /V)+ 

as (1 - < |x| < r^. 

Proposition 1.2.3. Let M = 1, = 0, /V x Z^ and B>Z\. 

(i) Then 

(1.2.12) Wj^ < Z\x\-^ 
and 

(1.2.13) pY < CZ'2\x\-'2 + CBZ^x\-'2; 
(a) There exist and r'^, 

(1.2.14) r^xe-izi, r;xfi-iZ„, 

stic/i i/iai l/l/g"^ B as \x\ ^ r^, l/l/g"^ ^ — |x| ^ r'^ and then pY = 
iffx> rrr,; 

(Hi) (1.2.12)-(1.2.13) become equivalencies (xj as \x\ < (1 — e)rm; 
(iv) 1/ X (Z - N)+bIz-1 < B and 

(1.2.15) 1/1/"^^ + 1/ X B\r^ - \x\y + r-\Z - N)+{r^ - \x\) 

- W^' X B\r^ - |x|)3 + r-2(Z - /V)+ 

as (1 - e)r„ < |x| < r^. 

Proofs of propositions 1.2.2 and 1.2.3. Proofs easily follow from equation 
and "boundary conditions" satisfied by w{r) where r — |x|: 



(1.2.16) w" + 2r'^w = P's{w + v), 

(1.2.17) w = r-'^Zrr, + 0{l) OS r^O, 

(1.2.18) w{r^) = -V, W{r) = vr'^ 

where u — — (Z^ — N)^?~^. □ 
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Corollary 1.2.4. Let M = 1, = and N x Z^. Then 

(V if \x\ > r'^ where r'^ x e^^Z^ as B > zl and r'^ x as 

B < cZi; 

4 

(a) As B < Zm the main contribution to both the charge and the Thomas- 

_ 1 

Fermi energy is delivered by zone {x : |x| x r^} with = Z^^ ; in particular, 

7 

then S^^ X S'^^ x Z^; further, in this case Wj^ x W'^^ in the zone 

4 

(^zuj Further, S^^ ~ as 6 ^ Z^; furthermore, in this case W^^ ~ 1/1/"'''^ 
in i/ie zone {x : |x| ^ rm}; 

4 

('zfj On t/ie ot/ier hand, as B > Zi , the main contributions to the total charge 

1 5 

and energy are delivered by {x : |x| x rm} and in particular Zm x BZmfm 
and 

(1.2.19) ^J^x ezlrix eizi; 

Recall that x as 6 < Z^^ and x B^sZ^^ as 6 > Z^. Note 
that proposition 24.3.5 of [Ivrll] (comparing W^^ for molecule with the 
sum of those for single atoms) still holds. Therefore we conclude that 

Corollary 1.2.5. (i) Assume that Z^ ^ Z for all m = 1, ... , M. Then all 
statements of corollary 1.2.4 remain true for M >2 with \x\ and Z^ replaced 
by £{x) and Z and r^, r'^, by r, r' , r* respectively; 

(a) In the general case global statements remain true, pointwise statements 
remain true without modification only as i{x) = im{x) ■= \x — Yml with 
Z^^Z. 

Also holds proposition 24.3.12 of [Ivrll] as it uses only super-additivity 
of t{p) and Te(p) is also super-additive (this follows from convexity of tb{p) 
and equality re(0) = 0). 

However there is a significant difference: if there is no magnetic field 
atoms really repulse one another on any distances and we can attribute it 
to either excessive positive charge as A/ < Z or their infinite spatial size as 
N = Z. However with magnetic field atoms have a finite size even as A/ = Z 
and they do not repulse one another on the large distances. In particular, 
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Proposition 1.2.6. Let N = Z and 

(1.2.20) \ym-ym'\>~rm + ?m' Vffl : l<m<m'<M. 
Then 

(1.2.21) £r{Z,y,B,Z)= {Zn^.Y^.. B, Z^) 

l<m<M 

and 

(1.2.22) p7(x,Z,y,e,Z)= ^ p7(x, Z^, y„, 6, ZJ. 

l<m<M 

As far as we know Theorem 1 of R. Benguria [Be] (see Theorem 24.3.8 
of [Ivrll]) has not been proven in the case of magnetic field; while it is not 
important for our arguments we still formulate 

Problem 1.2.7. Prove (or disprove) Theorem 24-3.8 of [Ivrll] in the case 
of the constant magnetic field. 



2 Applying semiclassical methods: M = 1 
2.1 Heuristics 

Let us consider first a mock proof of our main results; we deal here as if 
W^^ was very smooth which it is not the case; however later we will show 
that its smoothness is sufficient to employ arguments of Chapter 18 rather 
than those of Chapter 13 of [Ivrll]. We also will deal as if non-degeneracy 
conditions were satisfied leaving them also to more rigorous arguments 
below. 

2.1.1 Total charge. 

Consider 

(2.1.1) j e{x,x,u)ij{x)dx 

first with 7-admissible ip{x), where 7 < ei. 
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General airguments. 

The main part of semiclassical expression for (2.1.1) is of magnitude h'~^ + 
H'h'-^ X (^7^ + 6(73 with h! = 1/(C7) and fx' = B-f/C- Really, let us 
rescale x 1— )■ x/7 and r 1— )■ r/^^ which leads to h — 1 i-^ H' and 1-^ //'. In 
particular, as 7 x £ we get 

(2.1.2) C^^^ + BCf. 

Meanwhile the remainder in the semiclassical expression for (2.1.1) does 
not exceed Ch'~^+CijL'h'~^ x ^^7^+67^ (gaining factor h' in the comparison 
to the main part; here we need smoothness and as /i' > h'^^^ we need non- 
degeneracy) ; as 7 X ^ we get 

(2.1.3) c^e^ + Bf. 

Sure, we ignored the fact that H' < 1 does not necessarily hold even 
if 7 X £ but we believe that the contributions to the main part and the 
remainder of these zones will be less than of zone where this inequality holds, 
provided B <^ Z^. 

Finally let us sum with respect to ^-partition. 

Moderate magnetic field. 

Consider B < Zl first. Then for £ < Z'l we plug C = Z\i-\ into (2.1.2) 
and (2.1.3) resulting in 

(2.1.4) 0^^ zi£i + SZ5£i and Zi^B^ 

in the main part and in the remainder respectively and the summation 

4 1 

over zone {i < Z"?} results in the same expressions as £ = Z~3 i. e. 

1 2 2 2 

Z + BZ~3 X Z and Zs + BZ 3 x Zs respectively. 

Further, for i > Z~3 we plug ( — £~'^ into (2.1.2) and (2.1.3) resulting 

in 

(2.1.5) 01 r^ + Bi and r^ + Bf; 

11 3 

then summation over zone {Z~3 < £ < ? = 6" 4} results in Z + 64 x Z 

212 
and Zs+SsxZs respectively. 
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Strong magnetic field. 

4 1 

Consider B > Zs now. Then threshold Z s disappears and we sum ex- 

12 9 3 

pressions (2.1.4)q over i <? := Z^B^^ resulting in Z^B^^ ^- Z ^ Z and 

_,6 2 _2 _,1 ' -.2 _,1 

Therefore, 

(2.1.6) The total charge is min(A/, Z) (due to the choice of u) with the 
remainder estimate 0(max(Z3, ZsSs)) which is 0{Z3^ as S < Z? and 

o(ziel) as Z5 < e < Z\ 

2 1 

Remark 2.1.1. Remainder is less than the main part as Z^B^ < Z i.e. 
B < Z^. It means exactly that C,^>la,s^ = ? (as B > Zs). The same is 
true for all other semiclassical asymptotics below. 

2.1.2 Semiclassical D-term. 

Consider now the semiclassical D-term 

(2.1.7) CD{e{x,x,u) - pY{x), e{x,x,u) - pY{x)). 
General arguments. 

We do not have appropriate asymptotics for e(x, x, z/) in the case of magnetic 
field^^ but we apply Fefferman-de Llave decomposition (16.3.1 of [Ivrll]): 

(2.1.8) |x - y\-'{x, y) := |x - y\ip{^-^\x - y\) = 

where if E ^°°([1,2]). 

Therefore contribution of 6(z,7) x B[z','j) with 37 < |z — z'\ < 47, 
7 < ei{z) to such term does not exceed (C^7^ + 67^)^7~-^. There are x ^-^7"^ 
of such pairs with i{x) x i and their total contribution does not exceed 
C(C' + B)H\ 

Now we need to sum over 7"-^ c/7 which does not look good as leads to 
the logarithmic divergency but there is a simple remedy: we treat this way 

Unless we really assume that W is smooth and apply results sections 16.5-16.8 of 
[Ivrll]. 
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only pairs ti < \z — z'\ < i and apply for pairs with |z — z'| < ti pointwise 
asymptotics; then we get 

(2.1.9) C{e + Bfi'{l + {\ogBi/C)+); 

to get rid off this logarithmic factor we apply more delicate arguments 
similar to those of Subsection 16.9.3 of [Ivrll]. 

So, ignoring this logarithmic factor we conclude that the contribution 
of all pairs (z, z') with i{z) x i{z') x i does not exceed C{(^ + B^i^ while 
contribution of all pairs (z, z') with i{z) x £i 5^ £(z') x £2 does not exceed 

C{a + B){Ci + B)ilil{i, + i2)-\ 
Finally let us sum. 

Moderate magnetic field. 

4 11 

Consider B < Zs. Then summation over zone {ii 

5 1 

results in CZ3 and the same is true for summation over zone {Z^s < ii < 
B-\,Z I < ^2 < B--^}. 

Obviously, in such estimates, if there is a fixed number of zones, we do 
not need to sum over "mixed" pairs . 

Strong magnetic field. 

4 12 12 

Consider B > Zs. Then summation over zone {£1 < Z^B^^,i2 < Z^B^^} 
3 4 

results in CZ^Bs. 
Therefore 

(2.1.10) Term (2.1.7) does not exceed C max(Zi, Zist) which is CZi as 
B<Zl and CZIbI as zI < B < Z^. 

2.1.3 |Aa/ — i^l and another D-term. 

Consider two other non-trace terms in the upper estimate. 

Moderate magnetic field. 

As 6 < Zs we established the remainder in the total charge 0(Z3). Then 
using our standard arguments we conclude easily that |A/v — z/| = 0(Z) and 
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then 

(2.1.11) |Aa/-i^|- |N(z/)-A/| < CZi 
and 

(2.1.12) D(P^(l/l/J^(x) + Aa/) - P'B{Wr{x) + u), 

Pb{W7{^) + A/v) - P'siWrix) + ^)) < CZl; 
combining with the estimate of the previous subsubsection we conclude that 

(2.1.13) D(p,, - pY, - pY) <CQ = 0(Zi) 
as in (24.4.55) of [Ivrll]. 

Strong magnetic field. 

4 o 

Let now 73 < e < Then we estabhshed the remainder in the total 

2 1 

charge 0{Z^B^) and for the semiclassical D-term we established estimate 

3 4 

0{Z^B^) and therefore to estimate 

(2.1.14) \Xn - N\ < CzIbI 
as well we want to prove that 

(2.1.15) |Aa/ - z/| = o(z^et). 

Note that |z/| < Z? x Z^B^ < CB. Therefore if |Aa/ — < 2 1^^! 
conclude that 

(2.1.16) I j (P^(l/1/J'(x) + Aa,) - P'B{Wr{x) + i^)) dx\ > 

€\\n-u\B J {W + dx 

with integral taken over zone {x : W{x) + u > |Aa/ — z/|}. 

One can see easily that as|AA/ — z^|<e|z/| the right-hand expression of 

12 2 1 

(2.1.16) is larger than e|AA/ — ■ Z^B s and it must be less than CZsB^: 

12 2 1 

\Xi\i — v\Z^B^^ < CZ^Bs which implies (2.1.15). 
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Let us estimate the left-hand expression of (2.1.12). For this however 
estimate (2.1.15) is insufficient. We consider here only atomic case. Then 
using (1.2.15) one can prove easily that the right-hand expression of (2.1.16) 
is of magnitude 

|A/v — z^l ■ B?^ X {\u\?^^)^^ B^^ X |Aa/ — ■ B~^ 
provided |A — A/| < eu, where the selected factor is just integral of 

_ -- 2 1 

[B^z^ + \h'\r^'^)^^ dz (appearing due to (1.2.15)). Comparing with ZsB^ we 
conclude that 

(a) If \u\ > QZ-^oB^o (= Ziei x Z-^Sk) then 

(2.1.17) |Aa/ -iy\< C\u\'iZ-T5Bi> 

which is less than e|z/| and coincides with (2.1.15) as (Z — A/)+ x Z; 

(b) If |z/| > CiZ^^B^ then |Aa/ — < C2Z^^ B^ . 

In the former case one can prove easily that the left-hand expression of 

3 4 

(2.1.12) does not exceed CZsB^. 

In the latter case (as in Subsection 24.4.2) of [Ivrll] we consider Thomas- 
Fermi theory with z/ = i.e. N = Z and also prove that that 

3 4 

(2.1.18) The left-hand expression of (2.1.15) does not exceed Q = CZ^Bs. 

13 2 

In particular we slightly improve estimate (2.1.14) to \u\3Z5B3 as well 
(if (Z-A/)<Z). 

Therefore in our framework we estimated all non-trace terms in the 

3 4 

upper estimate by CZsSs and therefore "proved" estimate 

(2.1.19) D(pvu - pY, pw - p7) <cq = o(ziet). 

2.1.4 Trace 

Consider now Tr^ [Ha w — u). This term is of magnitude J {(^ + B(^) dx and 
one can see easily that it is x Z^ as 6 < Zt and x siZs as Zt < 6 < Z^. 
Meanwhile consider the remainder. Again for simplicity consider only 

1 5 

atomic case. As 6 < Z contribution of zone {x : i{x) < Z~i} is 0(Z3) (we 
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need to include Scott correction term in the main part) while contribution 
of zone {x : i{x) > Z~3} does not exceed 

(2.1.20) C J {C^ + BQr^dx 

5 

taken over this zone and it is x Zs as well. 

As Z < 6 < S'^ contribution of zone {x : i{x) < b .= S^sZa} is 
0{b~2Z2) = 0{Z3B3) and we need to include Scott correction term. Mean- 
while contribution of zone {x : £{x) > b} does not exceed integral (2.1.20) 

4 1 3 4 

taken over this zone which is :< + Z^B^ where the last term coincides 

4 5 4 

with estimate for (2.1.7) as B > Zs and docs not exceed CZs as B < Z^. 

Finally, as Z^ < B < B^ we need to reset b = Z~^ as h — 1/(C^) becomes 
< 1 inside; then we do not need Scott correction term and contribution 
of zone {x : £{x) < b} to both main part and remainder do not exceed 
C /(C^ + BC^) c/x X Z^ + e X e and contribution of zone {x : i{x) > b} to 
the remainder docs not exceed integral (2.1.20) taken over this zone which 

3 4 

results in CB + CZsBs and the second term dominates due to assumption 
B <^ Z^. Thus we arrive to 

(2.1.21) The main therm in Tr~{HA,w — ^) is of magnitude Zs as 6 < Z? 

2 9 4--) 5 

and B^Zs as Zs < 6 < Z while remainder estimate is 0{Z3) as B < Z, 

0{zIbI) as z < b<zI, 0{zIbI + zIbI) as Z5 < e < Z^. 

As 6 < Z* we need to include into main part Scott correction term. 

2.1.5 Discussion 

Now let us formulate our expectations 
Remark 2.1.2. We expect 

(i) Estimate (2.1.13) as B < Zl and (2.1.19) as Zi < 6 < Z^; 

4 

(ii) Furthermore as for B < Z^ the main contribution to all terms needed 
to derive this estimate is delivered by zone {x : £(x) fti Z~3 j and effective 
magnetic field is ji = Bi/( ^ BZ~^ we expect improved to "o" (or better) 
estimate (2.1.13) as 6 < Z and a > Z^i 

Recall that a = mini<m<m'<M IVm — Ym'l is the minimal distance between nuclei. 
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(iii) Statement (ii) should be also true for the trace term; however then we 
need to include Schwinger term; 

(iv) The remainder estimate for the ground state energy is maximum of 
the remainder estimate for non-trace and trace terms; therefore we expect 
the same remainder estimate as in (2.1.21); statement (ii) should be also 
correct for the ground state energy; however then we need to include both 
Schwinger and Dirac terms; 

(v) We expect the described remainder estimate of the trace term and the 
ground state energy if a is large enough; otherwise it should contain term 
0(a~2Z2) as 6 < Z4 and a > (and in this case we include Scott 
correction term). 

4 4 

Remark 2.1.3. The other difference between cases B < Zs and B > Zs is 
that /ih = BC,~^ 1 in the former case as £(x) < ?; however in the latter case 
it happens only as ^{x) < B~^Z but in the zone 
an opposite inequality holds. 



2.2 Smooth approximation 

An approach described in Subsection 2.1 hits two obstacles: non-smoothness 
of W^^ and its possible degeneration i.e. VWj^ is not disjoint from 0. 
However non-smoothness of W^^ is due to the non-smoothness of Pb- So 
we want to consider first the zone where we can just replace Pb( 1/1/ + v) by 
PiyV + u) and therefore by some smooth function 1/1/ which does not 
necessary coincides with W^^ . 



2.2.1 Trivial arguments 

Obviously we can do this as an effective magnetic field — Bi/( < 1. In 
this case we do not need assumption W + u >^ C.^ and therefore we can 
take C = as 6 < Zs and i > Z~3 and ( = Z2£~2 in all other cases. 
Therefore zone in question is 



(2.2.1) Xi := {x : £(x) < 

with ?i 



as 1 < e < Z, 

e-iZ5 as Z<B<Z\ 



22 



In this zone Xi for such modified W we can unleash the fuU power of 
the same smooth theory as in section 24.4 of [Ivrll] and prove the following 

Proposition 2.2.1. Let 1 < B < Z^. Then 

(i) Contribution of zone Xi defined by (2.2.1) to 

(2.2.2) j (^e(x, X, i^) - P'{W{x) + z/)) dx 

2 

does not exceed CZ^, its contribution to 

(2.2.3) D(^e{x,x, u) - P'{W[x) + u), e(x, x, u) - P\W{x) + z/)) 

5 

does not exceed CZi and its contribution to 

(2.2.4) j [ei{x, X, v) + P{ W{x) + z/)) dx - Scott 

does not exceed CZ^ + Ca^^zt + CB^Z^ ^\ 

(a) Further, as B <^ Z and a ^ Z^5 we can recover for these contributions 

2 5 5 

estimates CZ^, CZ^v and CZiv with 

(2.2.5) v:= Z-^ + [aZ'^)-^ + [BZ-^Y 

respectively where expression (2.2.4) should be modified to 

(2.2.6) j (ei(x, X, v) - P{ W{x) + z/)) dx - Scott - Schwinger; 

Furthermore in this case contribution of Xi to 

(2.2.7) \ j ^' ^' ^^^^ ~ 

does not exceed CZ^~ . 



As a < Z -"^ we skip Scott and reset a — Z in the remainder estimate which 
become CZ^. 
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Remark 2.2.2. (i) So far we should use P(.) instead of Pb{ ) but we will 
prove that the same results would hold for Pg as well; 

(ii) In the next Subsubsections we expand this zone to one defined by 
fji < h~3 but for trace term we still need a separate analysis as < 1; 

(iii) The same estimates hold as we replace in all expressions (2.2.2)-(2.2.6) 
P by Pg; 

(iv) We assumed that B < as otherwise h> 1 not only in Xi but even 
in {x : > B}; 

_ 1 

(v) Note that as ?i> [Z — N)_^_^ this zone (and the whole analysis) could 
be cut short as outside zone in question W + u > 0. From paper [IvrlO] we 
already know how to deal with such irregularities; 

(vi) We need to assume that a > Z~5 and to include the second term 
(aZs)"*^ in the definition of v only as we estimate the trace term (2.2.4). 

2.2.2 Formal expansion 

Now we want to expand zone Xi. Note first that 

(2.2.8) P^(l/V + i/) - P'(M/ + i/) = 0(ei) 
and 

(2.2.9) Pb{W + u) - P{W + u) - KiB^{W + = 0(6i). 

Really, one can consider P'^iw) and Pb{w) as Riemann sums for integrals 
P'{w) and P{w) respectively; see Appendix 8.C for details. 

However under non-degeneracy assumption |V1/1/| x (^'^£~^ we can do 
better with the integrated expressions. 

Proposition 2.2.3. Assume that in B{z, 7) 

(2.2.10) |V"H/| < CaCV"' Va : \a\ < K, 

(2.2.11) |Vl/l/|>eCV^ 

-^"^^ Or even to /x < h~i under non-degeneracy assumption (2.2.11) with — £, in 
particular, in the atomic case. 



24 



and 

(2.2.12) B < C^. 

Then 

(2.2.13) J (/)(x) (P^(l/l/(x) + ~P's{W{x) + u)) dx = 0{B^CWl 

(2.2.14) j 0(x)(Pe(l/l/(x) + P^(l/l/(x) + i/)) dx = 0{B\-^^^) 
and 

(2.2.15) D (0(x) (P^( H/(x) + i/) - P^( H/(x) + v)) , 

0(x)(p^(i/i/(x) + z/) - p^(i/i/(x) + v))) = o(e\-V) 

(2.2.16) Pb{w) := P(i/i/) + {kiP"{w)B^ + ^aP'^fi^) ■ (l - (^(i/i^/S)) 
w/iere e <^°°([-2, 2]), = 1 on [-1, 1]. 

Proof. Rescaling x i-> X7~-^, w i-> m/C~^ and therefore B /3 = B(~^ 
one can reduce the case to 7 = C = 1, /3 < 1. Then estimates (2.2.13) 
and (2.2.14) are trivially proven by (multiple) integration by parts which 
integrates Pp on each step increasing its smoothness^^^ 

To prove estimate (2.2.15) we apply decomposition (2.1.8). Integration 
by parts shows that (2.2.13) with t-admissible function cj) is 0{/3^t2) as t > /3 
and therefore contribution of zone {(x,y) : |x — y| x t} is 0{P^t^ x t~^) 
and therefore the total contribution of zone {(x,y) : |x — y| > /3} is 
0(/3^). Meanwhile a total contribution of the zone {(x,y) : |x — y| < /?} is 
0{f3^ X 0^). □ 

Therefore we expect that the zone Xi defined by /x < 1 could be expanded 
to the zone defined by yU < Z?^? 1°) or even larger^^-'; furthermore, under 
assumption |V1/1/| x we can define X[ by < /)^5 or even larger^^\ 

In fact one can prove then estimates 0(/3^) but adding correction terms ^ HkP'^''- 
However this improvement is not carried on to (2.2.15) in full. 

Wc do not need for each £ have a sharp remainder estimates but need only them to 
sum to a sharp estimate. 
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2.2.3 Expansion: justification 

Now however we need to deal with e(x, x, u) rather than Pg( l/l/(x) + u) (etc). 

Proposition 2.2.4. Assume that in 6(z,7) conditions (2.2.10), C7 ^ 1 

and 

(2.2.17) B < cC^(C7)"^ 
are fulfilled. Then for •j-admissible (f) 

(2.2.18) j </)(x)(e(x,x,z/)-P^(l/l/(x) + z/))dx=0(CV), 

(2.2.19) j 0(x)(ei(x,x,r)-Ps(l/l/(x) + z/))c/x=O(C'7) 
and 

(2.2.20) D [(j){x) (e(x, x, u) - P^( W[x) + z/)) , 

0(x)(e(x,x, z.) - P's[W[x) + z.))) = 0(CV). 

Proo/. Estimates (2.2.18) and (2.2.19) are due to Chapter 13 of [Ivrll]. 
Really, rescale x i— )■ x'-y^'^, r r^^^ and /j = 1 = 7~"'"C~''", B fi = 

To prove (2.2.20) let us apply decomposition (2.1.8); then according 
to (2.2.18) J{t) (which is expression (2.2.18) with t7-admissible 0t) does 
not exceed C(^^7^t^ as long as tC7 > 1; therefore contribution of zone 
{(x,y) : |x — y| X t} into the left-hand expression of (2.2.20) does not 
exceed C{C^h^f x ^^7-1 x CCV- 

Then summation over t > = B^-^'y~^( returns C(^^^ j t^^ dt x 
C(^'j^ log/i (we assume that /i > 2; case yU < 2 has covered). So, the total 
contribution of zone {(x,y) : |x — y| > fi"^} does not exceed C(^j^ log/x. 

Let us get rid off logarithmic factor. Returning back to 6(z, t) stretched 
to 6(0, 1) one can see easily that conditions of proposition 13.6.25 of 
[Ivrll] are fulfilled as weU with T = m\n[t-\ h~h^) and thus \J{t)\ < 
C{ht-^y^T~^ < Ch-^t^{t^ + hh-^). Plugging into (2.1.8) we get 

ch-'r^ C t-Ht' + hh-') dt X ch-'r^ = ccV- 
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On the other hand, in zone t < ji we use the trivial estimate 

e(x, X, u) - P'{ W{x) + u) = 0(/iC V) 

(due to simple rescaling x h> fix) and its contribution to the left-hand 
expression of (2.2.20) does not exceed C(/iC^7^)^ x Ai~^7 ^ (^C^l^- D 

Combining with estimates (2.2.8) and (2.2.9) we arrive to the statement 
(i) below; combining with Proposition 2.1.7 to the statement (ii): 

Corollary 2.2.5. Assume that in B{z,'~f) conditions (2.2.10) and Q'-y > 1 
are fulfilled. Let be -admissible function. 

(i) Let 

(2.2.21) B<cC^-l; 
then 

(2^2^22) /,(x)(e(x,x,.)-P',(M.M + .)).x^O(CV), 

(2.2.23) J 0(x) (ei(x, x, r) - ~Pb{W{x) + z/)) dx = 0(C'7) 
and 

(2.2.24) D (^0(x) (e(x, x, u) - P'^i W{x) + z/)) , 

0(x)(e(x,x, u) - P's{W{x) + u))) = 0(CV); 

(ii) Let assumption (2.2.11) be fulfilled and 

(2.2.25) B<cC'j-l. 
Then (2.2.22) -(2.2.24) hold. 

2.3 Rough approximation 

Unless our analysis has been cut short with ?! > (Z — A/)^_^, we need to 
consider zone {x : £(x) with redefined ?i, so that this zone is described 

1 3 

by /i > /i^3 or fi > in the general or non-degenerate (i.e. satisfying 
assumption (2.2.11)) cases respectively. 

In this zone we consider ^^-mollification with e <^ 1. In contrast to 
potentials considered in Chapter 18 of [Ivrll] function Wj^ is more regular. 
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2.3.1 Properties of mollification 

First, recall regularity properties of Wj^: 
Proposition 2.3.1. Wj^ have the following properties: 

(2.3.1) |V"H/J''(x)| < c«CW^£(x)-l"l Va : \a\ < 2, 

(2.3.2) |V"(l/l/J''(x) - W^^{y)) \ < 

CoBi{x)-'2\x - + CoC{xfi{x)-^\x - y\ 

y\a\ =2 Vx,y : |x-y| < e£(x) 

where we recall 

(2.3.3) CM = min(zi£(x)-i,£(x)-2) as 6 < 

(2.3.4) C{x) = Zh{x)-'^ as B>ZI; 

Proof. This proof is rather obvious corollary of the Thomas-Fermi equation 
(1.1.3). See also arguments below. □ 

Let us consider B{z,£{z)) with ^ B and rescale x H> xi~^, W i-)- 
w = (~^{W + u) (where we included u for a convenience). After such 

5 

rescaling w E uniformly, but there is more: Thomas- Fermi equation 
(1.1.3) translates into 

(2.3.5) Aw = fP'^{w) = fCP'{w) + fCiP'^iw) - P'{w)) 

with /3 = B(~^; we used that Pb{^) is positively homogeneous of degree 3 
with respect to (1/1/, 6). 

Note that parameter rj := (i'^ < 1 and x 1 if and only if 6 < Zs and 
i > Z~3 (in which case ( x £~^). 

Also note that the first term and the second terms in the right-hand 

5 1 

expression of (2.3.5) belong to ^2 and ^77^2 respectively with the uniform 
classes and 

(2.3.6) = PCf = BC'f, V ■■= a' as /3 < 1 

9 5 

Because of this w E ® /Sr]^^'^ again with uniform classes. Iterating we 
conclude that w e © (5%]^"^ with arbitrarily large exponent K. 
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On the other hand, ii B > (i.e. /3 > 1) without invoking Pg one can 

5 

prove easily that w e 77^2 with 

(2.3.6) ' T] := PCi' = BC'f as /3 > 1 
Therefore 

(2.3.7) w e ® /3r]^^ with arbitrarily large exponent K as /3 < 1 and 
w e rf€^^ as ^ > 1 

and one can see easily that 

(2.3.8) Parameter t] = BCT^I^ is 0(1) and 77 x 1 iff either B < Z\ and 
i X or 6 > Zs and i x 6~5Z5 (i.e. near border of supp pg"^, uncut by 

Remark 2.3.2. It may seem strange to define 77 diflferently as /3 < 1 and 
^ > 1 but there is a good reason for this as we consider M > 2. Anyway, 77 
is the magnitude of the right-hand expression of (2.3.5). 

Proposition 2.3.3. (i) Let he a e -mollification ofw with e < min(/3, h^). 
Then as /3 < 1 

(2.3.9) |V"(m/ - M/e)| < Ca/377£i-l"l Va : \a\ < 2, 

(2.3.10) \Pp{w) - Pp{w,)\ < cprje'^ 
and 

(2.3.11) iP'piw) - P'p{w,)\ < cfi-q'^ei + c^77£i; 

(ii) On the other hand as 13 > 1 the right-hand expressions of (2.3.9) -(2.3. 11) 
should be replaced by similar expressions albeit without (3: 



(2.3.9)' 


\V{w- 


We)\ < Ca77£i-I"l 


Va : 


\a\ < 2, 


(2.3.10)' 


\PM- 




Va : 


\a\ < 2, 


and 










(2.3.11)' 






Va : 


\a\ < 2. 
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(in) Further, under assumption |Vi/i/| x 1 m both cases 

(2.3.12) I J (t>{x){P^{w) - Pfs{w,)) dx\ < cr]ei, 

(2.3.13) I J (t>{x){P'i3{w) - P'piwe)) dx\ < cr,s-^ 
and 

(2.3.14) D(0(P;,(m/) - P'p{w,)), <p{P'^{w) - P'^{w,))) < cr,hl. 

Proof. Proof of (i) is trivial; in particular we note that 7752 < /3. 

Proof of (iii) is also easy as then is different from w on the set of 
measure x /3~-^£ as /3 < Cq and on the set of measure x e as /3 > Q. 
Actually w is uniformly smooth as /3 > 1 and l[x) < e? and we do not need 
any mollification here. 

One definitely can improve estimates (2.3.12)-(2.3.14) but we do not 
need it. □ 

Consider now the analytical expression and estimate the semiclassical 
error. From now on until the end of the subsection we assume that M = 1 
to avoid possible degenerations. 

Remark 2.3.4- Note that we can reduce operator to a canonical form as 
e > C{fi^^h)^ \ \ogfi\ (see Section 18.7) of [Ivrll] but here we will have a 
much better estimate as we will take e > hi~^ . 

2.3.2 Charge term 

Let us consider the charge term i.e. expression j e(x, x, v) dx = (Tr 9(z/ — /-/)). 
Regular zone. 

Then Section 18.9 of [Ivrll] implies that as 

(2.3.15) l/l/ + z/xC^ 
and 

(2.3.16) |Vl/l/|xC^r^ 

contribution of the ball 6(x, £(x)) to expression (2.2.2) does not exceed 
C(l + fih)h~^ X C(^i^ + CBi^ exactly as in the mock proof. 
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Then summation with respect to ^-partition in this zone results in CBs 
as e < Z, CZi as Z < fi < Zt and Cfi^zi as Zt < 6 < Z^. 

Remark 2.3.5. (i) Condition (2.3.15) is fulfilled as £(x) < er. 

(ii) Further, as M = 1 both conditions (2.3.15) and (2.3.16) are fulfilled as 
|x| < (1 — e)?m (we pick up ym = 0). 

Border strip. 

Now we need to consider the contribution of the border strip 3^ := {x : 
l{x) < e} with 7(x) 



(2.3.17) r 
and 

(2.3.18) C 



and scaling we get 

Let us consider first the case u — 0. Then both conditions (2.3.15) and 
(2.3.16) are fulfilled albeit with £i = 7(x)r and ^^(x) = (^{x)^ instead of £ 
and (. 

Thus as <J7r > 1 (i.e. 7 > 7 := /i^) the contribution of the ball 
6(x, 7(x)r) to the remainder does not exceed Cjihr^'^'^ and therefore the 
total contribution of zone iVi := {x : 7 < 7(x) < er} to the remainder does 
not exceed 

(2.3.19) C^h-^ J -f{x)-^dx X Ciih-^l \ogh\ = CBr^\ \ogh\ 

which is 0(Zi) as long as 6 < Z5(logZ)~^. 

Further, the same approach works as < ('^^^ x (^h x (r~^ which is 
equivalent to (Z - N)+ < ( (then |VI/1/| x <;^£^^ as 7(x) > 7) and also if 



e(r — |x|)r ^ and r := r^. Here ^ x r, C x C with 
{Z-NyJ B<{Z-N)l, 



B-\ 
Z'^B' 




{Z-N)i<B<Zl 
zl<B< CZ^ 

B<{Z-N)i, 

(z-A/)f <e<z|, 
z| < e < cz^ 



31 



this condition is violated but \y\ < C^; in the latter case we need to pick up 

7 = 71 := kl^C"^- 

To get rid off the logarithmic factor let us consider propagation. Recall 

that it goes along magnetic lines i.e. as (xi,X2) remain constant. Let 
us consider propagation in the direction in which \x^\ increases (i.e. 7(x) 
decays); we do not need to consider zone iVi n {jxaj < Z^'^r} as contribution 
of this zone (2.3.19) is o{B?'^). 

One can see easily that we can follow dynamics which does not return 
for a time 7"j'(x) := 7"i(x)(7(x)/7)^ where 7"(x) x i^^i^^^ x ?C^^^^^ is a 
time required for this dynamics to pass though 6(x, ^i(x)). Therefore one 
can replace (2.3.19) by 

(2.3.20) Cnh-^ I 7S(x)-^-^ dx x C/x^-^ = CB?^. 

J{x:7(x)>7} 

Further, as > Br we need also to consider zone ■= {■'^ ■ 7(-'^) < 7}- 
In this zone we take ii = ii = and <^ = ^ = (1^^17)2 with ii<; > 1 and 
since |V1/I/| x ^'^(-i^, contribution of 6(x, £i(x)) to the remainder does not 
exceed CBl\ and the total contribution of 3^2 docs not exceed CB?^ which 
what exactly we got for zone after we got rid off logarithm. We take 
mollification parameter e = ^~^Z^ 

Furthermore, zone 3^3 = {x : |x| > r + 4} is classically forbidden. So we 
can take here 

(2.3.21) £i(x) = e(|x| - r), ^{x) = m\n{^I^hi{x)KW\'^) 

and prove easily that its contribution also does not exceed CB?^. 

Returning to the case < C we see that the contribution of zone 
3^2 to the remainder does not exceed CB?^ because effective semiclassical 
parameter here is /i* = 1 and non-degeneracy condition is of no concern for 
us. Wc take mollification parameter e = q~^Z^ 

Moreover, we can modify W in 3^2 (make it negative there) so that this 
zone would be classically forbidden with £i, <; defined by (2.3.21) with 
replaced by (. 

4 

Finally in the case B < {i. e. B < C{Z — N)^ we can apply the 
above arguments with 7=1 and arrive to the same result. Therefore we 
proved in all cases 

^•^^ One can sec easily that the resulting errors in the expressions (2. 2. 2) and (2.2.3)- 
(2.2.4) will not violate our claims. 
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(2.3.22) The total contribution of the border strip y to the remainder in is 
0{Br^) which does not exceed CB2 as 6 < Zs and CBsZs as Zs < S < Z^. 

Conclusion. 

As < B < we need to estimate also contribution of the inner core 
Xq := {x : £(x) < CZ~^}. By means of variational methods we will prove 
(see 8.B.2) 

(2.3.23) As Z^ < B < Z^ the contributions of Xq to both J e{x,x, u) dx and 
/ Pe(l^W + ^) do not exceed CBZ~^. 

Then we arrive to the following 
Proposition 2.3.6. Let M = 1. Then 

(i) For constructed above potential W expression (2.2.2) does not exceed 
Czl + CBiZl; 

(a) As B < Z expression (2.2.2) does not exceed C{B + l)'^Zi~'^. 
2.3.3 Trace term 

Let us consider the trace term i.e. expression J ei(x, x, u) dx = Tr((/-/ — 
Regulatr zone. 

Here again let us consider first zone where |x| < (1 — e)r. Then the 
contribution of B{x,i{x)) to the Tauberian remainder^^^ does not exceed 
C('^{h^^ + /u) X C(^£ + CB(i as in the mock proof and the summation over 

5 14 3 4 

zone results in CZs + CSsZs + CZsfis. 
Border strip. 

Again in zone contribution of B{x, 7(x)) does not exceed CS^^i and the 
summation over this zone returns 

(2.3.24) CB J ^£^2 dx 

We will consider a bit later transition from the Tauberian to magnetic Weyl 
expression. 
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and plugging £i = ?7 and ^ = (^7^ results in CB^ as B < and CB^Z^ 
otherwise. The analysis of zone 3^o if there 3^2 = is also easy. 

Consider zones 3^2 and 3^o- The same arguments as before imply that 
their contributions to the remainder do not exceed CBT^^i^^ which is what 
we got before. 

Justification: from Tauberian to magnetic Weyl expression. 

Case fih < Cq. We need to prove that with the announced error we can 
replace the Tauberian expression by magnetic Weyl one. Note that the 
canonical form of (~^Ha,w as described in Sections 13.2 and 18.7 of [Ivrll] 
is 

(2.3.25) n = no+ 

/i^^u;i(xi, fi~~^hDi, X3) + fj,~^uj2{xi, fi'^hDi, X3){x2 + ^^h^D^) 
+ ^'^hW3{xi, fl'^hDi, X3) + 0(/i"^n(7 + /i"^)"5 

with 

(2.3.26) no = h^Dl - (x| + i^h^Dl ± iih) + M/(xi, ii-^hDi, X3) 
and 

(2.3.27) 7 = e min — 2_/'yu/i| 

j 

where we used the fact that w e ^hf^^ + ^^^h = ■ ^h. Here we 
have signs "+" and "— " on q/2 of the diagonal elements equally. 
Then the Tauberian expression is 

(2.3.28) const ■ /i/i"^ j^i^ - 2i>^ - A^^V - 2jy.'^huj2)\^ 

{i) + + 2Jfx-^hip2) dx 

where term with J = enters with the weight | and an error does not exceed 

Ch-^i^fx-"^ + fx-^h j (7 + ^-^)-5 c/x) X Cfi-h-^ 

because an integral does not exceed C/x^/x/i)"-'^; since /i > h~5 this error 
does not exceed Ch~^ which is better than 0{h~^). 
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On the other hand, if we consider the difference between (2.3.28) and 
the same expression with wi = = V^i = "^2 = and consider it as a 
Riemannian sum and replace it by an integral we get Gfi~^h~^ with an error 
not exceeding Cfi~^{fih)2h~^ which is even less. Therefore (2.3.28) becomes 

3 

and comparing with the result as /i x /i^5 when we get the same answer 
albeit with G = we conclude that G must be 0. This concludes the 
justification in 

Case fih > Cq. In this case we need a simplified version of (2.3.25) "H = 
"Ho + 0{fi~^h) and we need to consider only j = and replacing "H by "Hq 
brings and error C^h^^ x ^^^h = 0{hr^). This takes care of X2 and after 
scaling of y. 

Conclusion. 

As < B < we need to estimate also contribution of Xq = {x : i{x) < 
CZ~^}. By means of variational methods we will prove (see 8.B.2) 

(2.3.29) As 72 < e < the contributions of Xq to both / ei(x,x, u) dx 
and / Pe(l/l/(x) + u) dx do not exceed CB. 

Then we arrive to the following 
Proposition 2.3.7. Let M = 1. Then 

(i) For constructed above potential W expression (2.2.4) does not exceed 

czl + cb'^zI + czIbI; 

(a) As B < Z expression (2.2.4) does not exceed C{B + lyZ^^^ (but one 
should subtract a Schwinger term from the trace). 

2.3.4 Semiclassical D-term: local theory 

Unfortunately, we do not have non-smooth theory (cf. Section 16.7) of 
[Ivrll] here but actually we almost do not need it as singularities are rather 
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rare. Let us introduce a scaling function (2.3.27) and consider 

(2.3.30) Ja(z) = j (t)zAx){e{x,x,r)-Pp{w{x) + r))dx 

with (pz.xix) = 0(A~"'"(x — z)) and A < 7(z). Scaling x ^ A~-^(x — z) we have 
II ^ fi' = Xfi and h' = \~^h. 

Then, according to Section 13.4 of [Ivrll] 

(2.3.31) \Jx[z)\ < Ch'-^{1 + fi'h') X CX^h-^{l + fih) 
as long as X > h. 

Really, transition from the Tauberian decomposition to magnetic Weyl 
one in this case is easy: skipping all perturbation terms 0(/i~^) in (2.3.25) 
and also setting ipi = i'l = results in error 0[n~^ ■ fih~^ + jj,~^h~^) in 
(2.3.28)-hke expression albeit with the power | rather than | and without 
integration: 

(2.3.32) const ■ ijh'^^{w - 2j^h - fi'^ui - 2Ji2~^huj2)lx 

[tp + + 2j^^'^hii>2) 

Let us apply this estimate (2.3.31) to the Fefferman-de Llave decompo- 
sition. 

Case /i < C^fr^. 

(i) Consider first a pair (z, z') such that |z' — z"\ < eo7(z'); then also 
|z' — z"\ < eo7(z") and we take A = e|z' — z"|. 

Then in the virtue of (2.3.31) the total contribution to D-term of all such 
pairs belonging to 6(z, 7(z)), and with |z' — z"\ x A does not exceed 

(2.3.33) C7^A-^ X A^^ x X^h-\l + iihf x C-i^h-\l + ^ikf 

where C7'^A~^ estimate the number of such pairs. A""*" the inverse distance 
between them, and CX^h~^(l + /i/i) is the right-hand expression of (2.3.31). 

Then summation over A G (/i~"'",7) results in C7^(l + ^hYhr'^\ log(/i7)|. 
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Further, summation over all balls 6(z, 7) C 6(0, 1) with j{z) x 7 results 
in C{fih)~^'-fh~^ \ log(/i7)| as there are x [^h)~^'~f^^ such balls due to non- 
degencracy assumption |Vi/i/| x 1. Summation over 7 e {ii~^,iih) results in 
Ch-'^\\og{^i^h)\. 

As A < /i^"*^ we can apply standard non-magnetic methods without Fefferman- 
de Llave decomposition. Coefficients are smooth after scaling as long as 
e > //-^ 

(ii) Consider disjoint pairs {z',z") with |z' — z"\ > max(7(z'), 7(z")). Here 
estimate (2.3.31) is not sufficient and it should be replaced by 

(2.3.34) |J^(z)| < CX^h-^{l + i^h) 

as 

(2.3.35) 7 > ^i-^ 

Really, the shift for time T with respect to ^3 is x T provided |Vx3i^| x 1 
and this shift is observable if 7" x 7 > h^~^. Similarly, in the canonical form 
the shift for time T with respect to yU"-^^,- is x ijT^T provided |Vx,iv| x 1 
and this shift is observable if n~^T x 7 > ii~^h}-~^ . In both cases shift with 
T e (72, Co) is observable under assumption (2.3.35) and therefore we can 
extend T x 72 to T x 1. 

Note that for e > h^~^ assumption (2.3.35) is fulffiled automatically. Then 
contribution of each disjoint pair to D-term does not exceed 

Ch-\l + iihf^{z'f^{z"f\z'-z"\-' 
and the total contribution does not exceed 

Ch-\i + ^lhf jj \z' - z"\-^ dz'dz" X Ch-\i + /ihf. 

(iii) To shed off logarithm in (i) we need a slightly better estimate than 
(2.3.31). The same arguments as in (ii) result in 

(2.3.36) \Jxiz)\ < CX^h-^{l + ixh) ■ (1 + X-y/h)-^. 
Really, we just advance from time TxAto TxA(l-|-A7/^)''. 
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Then the same factor is acquired by the right-hand expression of (2.3.33) 
and the summation with respect to A G (^17""^, 7) results in C^^{1 + iJ,hyh~^ 
but summation with respect to A G {fJ>~^, h'j'^) results 

Cj%i + fihYh-\i + 1 iog(/x%-^)|). 

Further, summation over all balls B{z, 7) C 6(0, 1) with j{z) x 7 results 
in C{nh)~^'y{l + fj,hYh~^(l + \ log(/i/i7~-'-)|) and, finally, summation over 
J ^ l-ih results in Ch~^{l + ^Kf-. 

Note that in all cases perturbation terms in (2.3.25) and (2.3.32) result 
in the error not exceeding the announced one. 

Case n > Coh~^. So far factor (1 + ^h) was for compatibility only. Now 
it is important. 

Exactly the same arguments work as /i > CqH"^ with a minor modifica- 
tions: 

(a) 7(x) now is defined by (2.3.27) with 7 = and its upper bound is 1 
rather than /ih; 

(b) Also the number of 7-balls is x 7"^ rather than x 7"^; 

(c) A now runs from to 7 in (i), (iii). 

(d) We need to estimate contribution of pairs (z', z") with \z' — z"\ < h. One 
can see easily that e(x, x, r) < fih^^ and therefore the total contribution of 
these pairs does not exceed Cjji^fr'^ JJ \z' — z"\~'^ dz'dz" x Cji'^h''^ x /i^ x 

Therefore we have the following 

Proposition 2.3.8. As |Vm/| x 1 and e > hi~^ in 6(0, 1) and (j) G 
<^~(6(0, 1) 

(2.3.37) D(0(e(x,x,T)-P^(iv(x) + r)),0(e(x,x,i/)-P^(iv(x) + T)) < 

C(l + lihffr'^. 

Remark 2.3.9. One can see easily that one can select e > such that 
expressions (2.3.12), (2.3.13) and (2.3.13) will be respectively 0{h^+^), 
0{h}+^) and 0{h^+^). 
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2.3.5 Semiclassical D-term: global theory 
Regular zone. 

The above results allow us to consider a total contribution of zone X2 into 
semiclassical D-term. As before let us consider ^-admissible partition of 
unity there and apply it to Fefferman-de Llave decomposition. Then the 
total contribution of the elements which are not disjoint does not exceed 

(2.3.38) J2^n'{^ + BC'YCC ^ / (C' + B')i' i-'di 

n ' ' 

where (l + 6^~^)^ and are (1 -|- jiK) and hr'^ respectively and is a 
scaling factor. 

Then as = Zi"^ an integral equals to the value of the selected 

1 4 

expression as £ reaches its maximum, i.e. as £ = Z~3 for B < Zs and 

12 4 o 5 3 4 

e = Z-^B-i for Zs < e < and wc arrive to CZ 3 and CZ 5 B 5 respectively. 

On the other hand, as (^^ = an integral equals to the value of the 
selected expression as I reaches its minimum, i.e. ^ — Z~3 and only in the 

4 5 

case B < Zi and we arrive to CZs again. 

Furthermore, the total contribution of the disjoint elements does not 
exceed 

(2.3.39) 1^" - ^pr'(i + ^c') (1 + Bgyicxa - 

n.P 

J I (£ + /)"'(C' + B)e{C + B)e^ l-'d££'-^d£'. 

Then as — Z£^^ and C'^ — Z£'~^ an integral equals to the value of the 
selected expression as both £ and £' reach their maxima, and we arrive to 

5 3 4 

CZs and CZiBs respectively. 

On the other hand, as = and ('^ = (wc do not need to consider 
mixed pair) an integral equals to the value of the selected expression as both 

5 

£ and £' reach their minima, and we arrive to CZs. 

Therefore we conclude (combining with proposition 8.B.2 as Z^ < 6 < 

Proposition 2.3.10. Let M = 1. Then 
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(i) The total contribution of the zone {x : £{x) < (1 — e)r} to the semiclas- 

5 3 4 4 4 o 

sieal D-term does not exceed CZi and CZs Bs as B < and Zs < 6 < 
respectively; 

(a) As B < Z this contribution does not exceed C{B + lYZ^~^. 
Border strip. 

Border strip 3^ = {x : (1 — e)r < i{x) < (1 + e)r} is more subtle. Here we 
need to use the same ii{x) — eo{? — \x\) partition as before. 

_ 1 

Remark 2.3.11. y is already covered by our arguments if r x (Z — N)j^^ . 

Close elements. Consider first contribution of elements which are not 
disjoint. It is given by the left-hand expression of (2.3.38) with C, replaced 
by ^i(x) — r7(x) and ^(x) = C'~i{xY respectively. However since the layer 
{x : 7(x) X 7} contains x 7"^ elements the right-hand expression should 
be replaced by 

j B^T^-f -f-^d-f X B^?^ 

since < B\ so wc arrive to 0(max(6t, Zist). 

Meanwhile for 3^2 we have 7(x) = 7 < 1 and ^(x) — ^ = with and 
its contribution does not exceed what we got for y^. 

Disjoint elements. Consider contribution of the disjoint elements. It is 
given by the left-hand expression of (2.3.39) with £, replaced by 7 and <; 
respectively. Note that \zn — Zp\~^ x f-'^'j-^j'-^ where we sum with 
respect to all pairs with 7n x 7 and 7p x 7'. Therefore the right-hand 
expression should be replaced by 

(2.3.40) y r^7-^d77'-M7' 

which leads to Cr^B^ \ log(r~-'^7)p which differs from what we got before by 
a logarithmic factor. To get rid off it we will use exactly the same trick as 
in paragraph "Border strip" proving proposition 2.3.6 because considering 
disjoint pairs we consider the same objects as there. Then instead of (2.3.40) 
we arrive to 

J fB^r'i~'f\ l~^d-fi-^di 
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which results in CT^B^. 

Meanwhile for 3^2 "we have 7(x) = 'y <? and <;(x) — ^ — ^7^ with and its 
contribution does not exceed what we got for yi. 

Conclusion. Finally, analysis in the outer zone is trivial. Therefore we 
arrive 

Proposition 2.3.12. Let M = 1. Then for constructed above potential W 

5 3 4 

(i) Expression (2.2.3) does not exceed CZs + CZ^Bs; 

(a) As B < Z expression (2.2.3) does not exceed C{B + lyzl'^. 



3 Applying semiclassical methods: M >2 

Let us consider now the molecular case (M > 2). The major problem is that 
non-degeneracy condition |Vm/| x 1 is not necessarily fulfilled. Therefore 
we need to find an alternative approach to the zone where ix> . Recall 
that it consists of three zones: zone X2 '■= {/i/j < Cq, VV^^ > eoC^} zone 
A'3 := {fxh > Co, Wj^ > eoC^}^^'', and the (most difficult) boundary strip 
y^{Wj'<eoe}. 

3.1 Analysis in zone A2 

3.1.1 Scaling functions in zone ^2 
Step 1. 

We will use the scaling method here; the good news is that Wj^ is sufficiently 
regular for a proper rescaling. Recall that after we rescale in x i— )■ £~-^(x — x), 
r 1-^ C^T in B{x, ^i) with i = i{x), ( — Z2£~^, the rescaled potential 
w = C '^Wj^ satisfies in 6(0, 2) equation 

(3.1.1) Aw = r]P'^{w) with rj = < 1, p = /iH = BC^ < 1 

Only as 6 < CiZ^; this zone disappears for CiZ^ < B < . 
Only as Zj < e < Z3; this zone disappears for 6 < Zs . 
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and therefore in 6(0, 1) 

(3.1.2) w = -— I \x - z\-\P'Jw{z))(f){z) dz ^ w' 

An J 

with G ^o~S(0, 1) and w' G <^o°°e(0, |). 
Let us introduce a function 

j_ 

(3.1.3) 7o(x) = (^min - 2jP\^ + |Vi/i/|^ + |VV|^ + |VVV^) 

5 

Note that we cannot replace w' by w in the last term because w 
only. 

Proposition 3.1.1. 7o(x) zs a scaling function i.e. \x — y\ < e7o(x) =^ 
7o(y) X 7o(^)- 

Proof, (i) If 1/1/ belonged to (and wc would put w instead of w' in the 
last term of (3.1.3)) then we would just prove that |V7o| < c. Here we 
should be more subtle. We need to prove that if 

(3.1.4) , 2 min |m/ - 2jP\ < 7o^ \Vw\ < 7^, 

J 

(3.1.4) 34 \V^w\ < -fl |VV| < 70 

at point X, then at point y the same inequalities hold with 70 replaced by 
C7o Definitely this is true for (3.1.4)4 since w' is smooth. 

(ii) Consider \V^w\. Consider \Aaw{y) - V"m/(x)| with A„ = V" - |5yA, 
a — Then due to (3.1.2) 

(3.1.5) \Aa,yw{y) - Aa,,w{x)\ = 

\ri J (A„,x|x - z|-^ - A^,y\y - z|-i)P^(i/i/(z))0(z) dz\ + e^-fl 

where the last term estimates |AQ,yi/i/'(y) — Aa^xw'[x)\ and we used (3.1.4)^. 
Integrals are understood in the sense of the principal value (vrai) and 
ei = ei(e) ^ +0 as e ^ +0. 



We need to prove a bit of reverse as well; see (iii). 
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Note that the integral in expression (3.1.5) does not change if we add to 
P^(m/) any constant with respect to z. 

Let us consider first this integral over {z : \x — z\ > 270}, as > \x — y\, 
and note that this integral does not exceed 

I f\ 1-41 I 

r)\ / |x — z| |x — y| X 

(|Vm/(x)| • |x-z| + |V^M/(x)| • |x-zp + |x-z|i)^ c/z| < Ceirj. 

Consider now integral over {z : |x — z| < 270} 

J A«,,|x - z|-^ • {P'f,{w{z)) - P'^{w{x))) dz\ 

and note that it also does not exceed eirj. The same arguments work for this 
integral with x replaced by y but still integrated over {z : |x — z| < 270} (one 
needs to remember that | Vi/i/(x) — Vi/i/(y)| = 0(70) and | V^i/i/(x)— V^i/i/(y)| = 

o{ih). 

Further, the same arguments work also for these expressions integrated over 
{z : |y — z| < 270} and we are left with 

\r} J uj{x,y,z){P'^{w{x)) - P'^{w{y))) dz\ 

integrated over {z : |x — z| x 70} and u x 7^^ and one can estimate it 
by ei?7 easily in the same way. Therefore since |Ai/i/(x) — Ai/i/(y)| does not 
exceed tii] we conclude that |V^m/(x) — V^w(y)| does not exceed 61(77 + 7q). 

However (3.1.2) implies that rj < cjq since P'j^iw) x 1 in and therefore 
\V^w{x) - V^w{y)\ < 6172 in e(x,7o). 

Combining with (3.1.4)2 we conclude that |Vw(x) — Viv(y)| < ei7Q in 
6(x,7o); next, combining with (3.1.4)^^ we conclude that |m/(x) — M/(y)| < ei7o 
in e(x,7o). 

(iii) Therefore (3.1.4)^^ ^ are fulfilled in y G B{x, ejo) with 70 replaced by 
7o(l + Cei). Further, if we redefine 70 as the minimal scale such that 
inequalities {3.1A)^_^ are fulfilled in x, then (3.1.4)^^ ^ fail in y e B{x, e^o) 
with 7o replaced by 70(1 — Cei). Therefore with appropriate e > we 
conclude that | < 7o(-'<')/7o(y) < 2. 
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Obviously 70 x 7ooid where 7ooid was defined by (3.1.3) and therefore the 
same conclusion also holds for 7ooid- D 

Now let us reintroduce the scaling function 

(3.1.3)* 7o(x) = e(m\n \w - 2jP\^ + |Vi/i/|^ + |V^i/i/|^ + \V^w'\^^Y' + 

CqH^ + Qr}^. 

Then 

(3.1.6) |x — y| < 27o(x) =^ 7o(y) ^ 7o(-^) and {3.1.4:)^_^ hold (with some 
constant factor in the right-hand expression). 

Consider 6(x, 70), 70 = 7o(^)i and scale again x 1-^ 7q""^(x — x), r i->- 7^V 
and respectively 1/1/ M- 1/1/1 = lo^{^ ~ 27/3)), h ^ hi = hj^^. 

If 7o X ^3 then ^1 x 1 and we are done. If 70 x 7^2 then go to Step 3. 

Step 2. 

So, let rji <^ 1. Let us introduce scaling function in 6(0, 1) obtained after 
the previous rescaling 

(3.1.7) 7i(x) = e(min |i/i/i + 2(7 - j)/?^-^!' + |Vi/i/i|3 + | V'l/i/il^) 
Then 

(3.1.8) ,_2 min |i/i/i + 2ij-j)P%'\ < CqjI |Vi/i/i| < QjI 

J 

(3.1.8) 3 |VVi| < Co7i- 

Proposition 3.1.2. (i) 7i(x) is a scaling function: |x — y| < 27i(x) =^ 
71 (y) >^ 71 (^); 

(^u^ If Tji < eo ^/ien 

(3.1.9) |V2m/i(x) - V^w[{x)\ < £271, K = {w' -pho'. 
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Proof. Proof is similar but simpler than one of Proposition 3.1.1; it is based 
on the rescaled version of (3.1.1)-(3.1.2): 

(3.1.10) Am/i = rjiPd^ilo + W). VI = Vlo^ < 1 
and therefore in 6(0, 1) 

(3.1.11) wi = -^J\x- z\~\P'f,{wi{z)fo + 2j/3)0(z) dz + w[. 

□ 

Now let us reintroduce scaling function 

(3.1.7)* 71 (x) = e(m\n + 2(7 - j)/37o"T + |Vi/i/i|' + | ' + Qhl 

Then 

(3.1.12) |x — y| < 27o(x) =^ 7o(y) ^ 7o{x) and (3.1.8);^_3 hold (with some 
constant factor in the right-hand expression). 

Let us consider x G 6(0, 1) (it is a new point), 6(x,7i), 71 = 7i(x), and 
scale again x 1— )■ 7f ""^(x — x), r Ii^t cind respectively i/i/i i—)- 1/1/2 = 'li^wi, 

_5 

h2 = hiji ^ . 

2 

If 71 X hi then h2 ^ 1 and we arc done. If |Vm/2| x 1 we are done as 
well. 

Step 3. 

So, consider the remaining case |V^m/2| >^ 1- Then we introduce scaling 
function 

1 J 

(3.1.13) 72(x) = e(min \w2 + 2{J -JWjq'j.^I + IV1/1/2I') ' + Qhl 

Let us consider x e 6(0, 1) (it is a new point), 6(x, 72), 72 = 72(^)1 and 
scale again x 1-^ "72 ^{x — x), r i->- 7^'^t and respectively W2 ^ w-^ = 7^^1/1/2, 

1 

If 72 X hi then /i3 X 1 and we are done. If |Vm/3| x 1 we are done as 
well. 
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step 4. 

Finally, introduce 



(3.1.14) 73(x) = emin \w3 + 2Q - j)Mo'li^l2^\ + Chi 



J 



3.1.2 Scaling arguments in zone A2: Semiclassical 
N-term. 

Now we apply scaling arguments using scaling functions constructed above. 

We revert our steps. While we call 71-3 relative scaling functions let 
us introduce absolute scaling functions aQ^x) = 7o(->^), Oii{x) = 7o(-'^)7i(-'^)5 
ot2{x) = 70(^)71 (^)72(^), and a-i{x) = 7o(x)7i(x)72(x)73(x) 

We need first 

Proposition 3.1.3. Consider 6(0, 1) and assume that in it 

(3.1.15) |Vm/|x^, 
and 

(3.1.16) \V^w\ < ce. 
Let 

(3.1.17) 7(x) ■=emm\w -2^Jihj\e-'^ + hl 

j 

and 

(3.1.18) e>{K)-i'\ K = he-'^. 

Let if e <^o°°([-eo, eo]). Then for a < ^ := 7(x), 

(3.1.19) I J 0«(x)(e^(x,x,T)-P;,,^(Kx) + r))cyx| < 

with 

(3.1.20) e^(x,y,T) := cp{fi'Dl{t,h)-')e{x, y,T) 



So far we ignore the first scaling x i— > (x — x)i ^. 
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and Weyl expression 

(3.1.21) P'i3,<p{w + r) = const ^/3(m/ + r - 2jiJ,h)^ip{w + r - 2j>/i) 

wzt/i t/ie standard constant where for each x only one term is present in this 
sum. Here we take large s > as < ^ya and s — otherwise. 

Proof. The proof is standard and based on the standard reduction to canon- 
ical form, standard estimates for U{x, y, t) a Schwartz kernel of e'^* ^ 

(3.1.22) \f=t^n-^rj XT{t)<t>a{x)U^{x,x,t)dx\ < Ciih-W 
as 7 X 1 and 

(3.1.23) |F,_,,vy" Ux){xT{t)-Xf{t))U^{x,x,t)dx\ < 

with T = e'-f^ where U^p is defined similarly to (3.1.20). 

Here obviously wc can skip in (3.1.19) all perturbation terms in the 
argument and in 0^ transformed. □ 

3 

Then plugging into (3.1.19) a = 7 (= 7), we have factor {h^y^^Y in the 
second term. 

There are two cases: 9 < /ih and 9 > /ih. 

In the former case taking the sum over 7-partition of 1-clcmcnt we 
estimate the same expressions with 0i instead of 0-y by their right-hand 
expressions integrated over 7"-^ dj which return C/ih"^. 

On the other hand, in the latter case let A = iJLh9~^. Taking the sum 
over 73-partition of A-element 4>x by the right-hand expressions which return 
C/Li/lj ^A^. In this case summation over A-partition return C9h~^. 

In both cases we arrive to 

(3.1.24) I j 0(x)(e^(x,x,T)-P;,,^(iv(x) + T))c/x| < C9h-^ + C^xh-\ 

Applying this estimate after a;2-scaling we conclude that the left-hand 
expression with (j) — 4>a2 (in the non-scaled settings) does not exceed 
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C6h~'^a2 + C^h~^a\. Here the first term is 0{h~^a^ and the summation 
over 1-eIement returns 0(h~^). 

Consider the second term Cfih^^a^ = Cyu/i^-'-7o7i7|. Then summation 
over a2-partition of ai-element returns 

C^h-^llll j 72 X 72"^ ^ C/in-So7i(l + I log72|) 

where 7/^ is a minimal value of 7/^ over 7/(_i-element. However in fact there 
will be no logarithmic factor because in virtue of equation (3.1.1) there is 
a positive eigenvalue of Hess 1/1/2 of the maximal size (cf. Section 5.1.1) of 
[Ivrll]. Therefore, in fact we have Cfih^^'-fQ'-fl. 

Now summation summation over ai-partition of ao-element returns 

C/irSo' I ^2 X ^-3 ^ c^rS2(i + 1 log^.i). 

Finally, summation over ao-partition of 1-element returns 

Cfih-' J 7o'(l + I log 71 1) X 7o"' dx ^ Cfih-'%\l + \ log 71 1). 

where 7^ is an absolute minimum of 7/c. However 'Jq > rj and 'ji > f] and 
therefore expression above does not exceed Cfih~^ri~^(l + | \ogr]\). 

Remark 3.1.4- Recall that we estimated only cut-off expression. To calculate 
the full expression we need to calculate also the contribution of zone > 
lih}. However this is easy. 

Really, instead of (p{h^ / {fih)) consider ip'{h^Dl/e) with ip' G ^°°([1,4]) 
and e/i/j < 9 < 1. Without any scaling one can prove easily that such 
modified expression (3.1.24) does not exceed C9h~^. We leave easy details 
to the reader. 

Therefore plugging 6 = 2"fih and taking a sum over n = 0, ... , [| log2 fih\\ 
we get required expressions. Also note that in such expressions we need to 
consider perturbed argument w + fi~^ui + Jfi^^hw2 (all other terms which are 
0(/i~'* + fi~2fi) could be skipped and also a perturbed function transformed. 

Remark 3.1.5. (i) However we need to do get rid off these perturbations 
for 9 < iih}~^ only. Really, for 9 > fih^~^ we need canonical form only to 
study propagation and calculations could be performed without it. But then 
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getting rid off perturbation is trivial provided perturbation does not exceed 
Cjih}-^^ which is the case if 

(3.1.25) n > h--^-\ 

(ii) Note that in the smooth approximation contributions of Xi is always less 
than CZi-^S C max(Zi, ZifitjZ-'^i or C max(zt, zistjZ-'^i + CfisZ^-^i 
respectively with an exception of first two and only in the case of the 
threshold value between > Z~3~'^2. However in this case 77 > Z~^^ and the 
errors of the smooth approximation approach in fact are less than CZ^~^^, 
CZ3~^'^ as well. Therefore there are in fact no exception. 

(iii) It is important to have e < jih* and with e = hi we need ji > 

which is due to (3.1.25). 

Therefore we conclude that in completely non-scaled settings with = 

(3.1.26) I J 0(x)(e(x,x,T)-P^(l/l/(x) + T))| < 

cCe^ + CBic Hi + \\oge^C\) 

where the first term is Ch~^ and the second term is Cfj,h~^rj'H^ + I 'og^D; 
recall that x x Bi(~-^ and rj x £^(. In comparison with the 

non-degenerate case IVI/I/J*^! x C,^£~^ we acquired the last term. 
Assume first that 

(3.1.27) Z;„xZ Vm=l, ...,M. 
Then 

(i) As 6 < Zt, £ < Z~i we have ( = Z^£~^ and the right-hand expression 

5 1 

of (3.1.26) returns CZ£ + CBl^Z"^ and the summation with respect to £ 

12 2 

results in its value as £ = Z~3 i.e. CZ3 -|- CBZ~3 with the dominating first 
term. 

(ii) As 6 < Zt, £ > Z~5 we have ( = £~'^ and the right-hand expression 

of (3.1.26) returns C£~'^ + CB£^. We need to sum as long as jih < 1 i.e. 

11 21 
Z~3 < £ < B~i and the summation returns CZ3 + CB2 with the dominating 

first term. 
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(iii) As Zs < 6 < Z^, i < B we have ( = Z2i \ and the right-hand 

1 5 

expression of (3.1.26) returns CZ£-|- CBZ~^i^. Then summation results in 
CZ^B-^ + CZ^B-'^Z <Zl 

If assumption (3.1.27) is not fulfilled we can estimate in the first term of 

the right-hand expression of (3.1.26) parameter ( from above by min(Z2£~2 , £-2^ 

1 _ 1 

and in the second term from below by Cm = rr]'\n[Zm£m^ , i^^) BS £ = £„, := 

— Ym I and repeat all arguments above. 

Therefore we arrive to the statement (i) below. Furthermore, as 6 < Z 

3 3 1 

zone X2 is contained in {x : £{x) > 6~ib > Z~io (really, n > h~3 in A2) 
and we arrive to the statement (ii) below. 

Proposition 3.1.6. (i) As B < Z^ contribution of zone X2 to expression 

(3.1.28) J (e(x,x,z/)-P^(l/l/(x) + z/))dx 

2 

does not exceed CZs; 

(ii) As B < Z contribution of zone X2 to expression (3.1.28) does not exceed 

czl-'. 

3.1.3 Scaling arguments in zone A2: Semiclassical 
D-term 

We need to estimate 

(3.1.29) D(0,(e(x,x, v) - P^(l/l/(x) + i.)), 0,(e(x,x, 0) - P^(l/l/(x))) 

where 0a (■'^) is an a-admissible function. Again we revert our steps. 

Consider 6(x, 0:3) and apply Fefferman-de Llave decomposition. Then in 
the framework of Proposition 3.1.3 contribution of pairs B{x, a) and B{y, a) 
with Za < \x — y\ < Aa docs not exceed the right-hand expression of (3.1.19) 
squared and multiplied by where Ca^^ estimates the number of the 
pairs and is the inverse distance. At this moment we discuss a cut-off 
version of (3.1.29) i.e. with e<^(., ., .) and P'^n^^i-)- So, we have 
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Then integrating this expression with respect to a da with a < 73 we 
arrive to C/i^/i~^ (73 + h^js). 
Therefore 

(3.1.30) A cut-off version of (3.1.29) with a = does not exceed C//^/i~^ (73+ 

The first term here C/j.^h'^'y^al does not exceed C/jL^h'^al (recall that 
aj = aj-i'jj) and the summation over aa-partition of 1-element returns 

Consider the second term C/x^fi~^/l273'^2! summation with respect 

to tts-partition of a2-clcmcnt returns Cjj?-hr'^a\h\ J 73"^ c/73 < Cjj?-hr^a2 

2 

(really, recall that according to (3.1.14) 73 > ^13) and then the summation 
over Q;2-partition of 1-element returns Cfi^h'^. 

Consider 6(x, 0.2) and apply Fefferman-de Llave decomposition. There 
are two kinds of pairs: 

(a) those with |x — y | > t{a-i{x) + az{y)) for all (x, y) and 

(b) those with |x — y| < min(Q;3(x), Q;3(y)) for all (x, y). 

The total contribution of the pairs of the second type (i.e. summation is 

taken over all pairs of «3-elements in 6(0, 1)) as we already know is 0{fi^h~'^). 
Meanwhile according to the analysis in the previous subsubscction Scaling 
arguments in zone X2: Semiclassical H-term a contribution of one pair of 
kind (a) does not exceed 

c(/.- + /,/r^,^^,r^,^^,^^-^7S«ix 

" V ' 

at X 

{h-^ + /x/i- Vo"'7rS2"S3"'"'73)«l X |x - y|-i 

" ' 

at y 

where each of two first factors is just an estimate of the integral (3.1.24) 
calculated over corresponding argument. If we take the first term in the 
first factor and sum over «3-partition of 1-element we get only the second 
factor multiplied by ^h^^ and then summation was done in the previous 
subsubscction. Similarly we can deal with the first term in the second factor. 
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On the other hand, if we take only second factors and sum over pairs of 
ccs-subelements of the same Q;2-element we get 

Then summation with respect to a2-pai'tition of 1-element returns Cn^hr^. 

Consider now 6(x, ai) and apply here Fefferman-de Llave decomposition. 
There are two kinds of pairs: 

(a) those with |x — y | > e{a2{x) + a2{y)) for all (x, y) and 

(b) those with |x — y| < min(Q!2(-'<'), ct2(y)) for all (x, y). 

According to the above analysis the total contribution of the pairs of the 
second type (i.e. the summation is taken over all pairs of Q;2-elements in 
6(0, 1)) as we already know is 0{ii^h^^). Meanwhile according to the analysis 
in the previous subsubsection Scaling arguments in zone X2: Semiclassical 
U-term a contribution of one pair of kind (a) does not exceed 

c (r^ + /xrSo"'7r'72"')«3 x (^"' + /^r^7o-^7r'72"')«3 xk - rl"' 

" V ' " ' 

at X at y 

and here again we can "forget" about the first terms in each factor. Then the 
summation with respect to pairs of a2-subelements of the same ai-element 
results in Cfj^^h^^jQ^j^^al x Cfj^^h^'^al where we avoid logarithmic factor 
in virtue of the same positive eigenvalue of Hess w. Summation with respect 
to tti-admissible partition of 1-elemcnt returns Cix^fT^. 

Consider now 6(x, qq) and apply here Fefferman-de Llave decomposition. 
Again there are two kinds of pairs and the total contributions of the pairs 
of the second kind we already calculated and contribution of the pairs of 
cci-subelements of the same 1-element does not exceed 

C/.2r^7o '(1 + I Iog7i|)'«o < Ci/h-\l + I \og7^\fal 
and the summation with respect to cco-partition of 1-element returns 

c^^^h-\l + \\og7]\f. 

Finally, consider B(x, 1) and apply here Fefferman-de Llave decomposi- 
tion. Again there are two kinds of pairs and the total contributions of the 
pairs of kind (b) we already estimated while the total contribution of the 
pairs of kind (a) docs not exceed Chr'^ -\- CiJ?'fi~^'r]~^{l -\- \ logr^l)^ where we 
recalled the forgotten terms. 
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Again, this is estimate for cut-ofF expression. Going to uncut expression 
we repeat the same trick as before but as we deal with D-term we need to 
consider "mixed" pairs when one "factor" comes with 9 and another with 9' 
but then contribution of such pair does not exceed C {i'h~^)2 [u' h~^)2 . Easy 
details are left to the reader. 

Therefore returning to the original scale we conclude that the contribution 
of £-layer to (3.1.29) does not exceed 

(3.1.31) CCf + CBHC\'^ + I log^^Cl)^ 

which is the right-hand expression of (3.1.26) squared and multiplied by 
due to scaling. In comparison with the non-degenerate case iVl/Vg"^! x C^£^^ 
we acquired the last term. 

Assume first that condition (3.1.27) is fulfilled. Then 

(i) As 6 < Zt, £ < Z~i we have ( = Z^£~^ and expression (2.2.18) returns 
CZ^i + CB^i2Z~2 and the summation with respect to i results in its value 
as £ = Z~3 i.e. CZs -|- CB^Z~^ with the dominating first term. 

(ii) As 6 < Z5, ^ > Z~5 we have C, = £~'^ and expression (2.2.18) returns 
C£~^ + CB^£^. We need to sum as long as u/? < 1 i.e. Z"3 < £ < 6"? and 

5 5 

the summation returns CZs -|- C6? with the dominating first term. 

(in) As Zl < B < Z^, £ < B'^Z we have C = Z\£-'\ and expression (2.2.18) 
returns CZ^£-|- CB'^Z~2£2. Then summation results in CZ^B~^ + CB2Z < 
ZiBi. 

Sure, we need to consider also mixed pairs of the layers and their 
contributions are 

c{cH^+CB£c'^{i+\ \og£^c\)) X (C"^"+cef C'~^(l+| logf 2^'|)) x(£+£')"' 

and the summation with respect to £ and £' returns the same expression as 
above. 

If assumption (3.1.27) is not fulfilled we use the same trick as in the 
previous subsubsection. Therefore we arrive to the statement (i) below. 
Applying the same arguments as in the proof of Proposition 3.1.6 we arrive 
to the statement (ii) below. 

Proposition 3.1.7. (i) As B < Z^ contribution of zone X2 x A2 to expres- 

5 3 4 

sion (3.1.29) does not exceed C max(Z3, Zsfis); 
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(a) As B < Z contribution of zone to expression (3.1.29) does not exceed 

czl-'. 

3.1.4 Semiclassical T-term: zone Xi extended 
What is Xi extended? 

First, let us analyze the precise extension in the framework of D- and N- 
terms. For these terms wc have errors C^iihYrj^^h^^ = CB^(~2i'^ and this 
expression squared and multiphed by i.e. CB^(^~^i^ respectively, and 
both sum to their values as i — r. 

(i) As 6 < Z5 and r > Z^5 we have CB^r^ and CB^r^ and we want them 

2 5 7 

to be less than CZ? and CZs respectively. As 6 < Ze the first requirement 
is more restrictive and we set CB^r^ = Z? i.e. r = B^^Z^s and then 
/i = Br^ — B~sZi5 > Z6 and h — B~iZis > Z~3 so n > hr^. 

(ii) As Zl < B < Z^ the second requirement is more restrictive and we 
set B^Z~2ri = Zs, i.e. r = B^^Zz? and a — B~3Z9 and h — Z~2r~2 = 

4 23 3 1 

B9Z~27 and iJ,> (better than h~3). 

But ah < rj iS r > B^Z~^ i.e. B < Zi and this is less than Zs so we 

5 11 5 5 

test /X and H in this case: fi = Z^ and h — Zm and [i > h~^. As B > Zn 
we have error {^Kf'hr^ r^^ = B^r^ and we want it to < Zs, so r = Z^B^^ 

1 3 2 4 1 

and fj, = Bw and h = SibZ^s and we test it as = Zs when /j, — B^o and 
h = B^5, so exponent — ^ fits. 

As D < Zs we consider B^r = B^Zs i.e. r = B 25Z25, ^ = Ssoz 25, 
h= 65oZ~25 and exponent — ^ fits as well. 

When we can use the same method for T-term? 

As far as semiclassical T-expression is concerned an error of such approach 
in the localized and scaled settings would be C{fih)2h^^ which is 0{h~^) 
only as fj,> h~i. One can extend it to > /i~5~^ using the same trick as in 
Remark 3.1.5 but we need to do better than this. 

Note first that in fact an error does not exceed C{iiKf'rj~2h~^ x CB^i^(~2 
in the localized scaled settings (we skip Ch~^ as the result here is already 
taken into account). The simple proof is left to the reader. This is translated 
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into CB^£'^(~2 into unsealed settings. Summation with respect to £ < r 
returns its value as i = r. 

So we get CB^r^ as 6 < Zs and r > Z~3. Let B < Z. Then we want 
CB^r^ < CZi i.e. r = B'sZs which is greater than Z~3 as 6 < Zt. Then 
fj, = Br^ = B~sZ > Zs and h = r = B~^Z^ > Z~k and ji > h~^. 

As Z < e < Z5 but A- > Z-5 stiU we want CB'^r^ < CB^Z^^ i.e. 
r = B~i^Zi^ and we want it to be > Z~3 i.e. B < Zs. Then u = Br^ — 

3 4 1 8 4 1 3 

B~5Z5 > Zb and h — B~^Z^ > B~3 and /j, > h~s. It is not as good as 

3 

/J, > h~T . 

Then we use the smooth canonical form. In the operator perturbation 
terms have factors /i^^, etc and we can use the standard approach to 
get rid off ji^'^ < fih, so we need to consider only 

However let before scaling the second derivative of W be of magnitude 
9; then after scaling it becomes of magnitude 9' — 70716* and then the 
perturbation is of magnitude ^//~^ but contribution of the error will be (after 
we compare the true Riemann sum and the corresponding integral and their 

difference Cn-^un-^{finf{9')--2 x a? < C9h-^^o^^^^al < CrVi'^o;? 
where we used that 9 < C7071. Then summation over ai-partition of ao 
element returns Ch~^aQ and the summation over 70-partition returns Ch~^ 
as desired. Therefore we covered zone Xi for T-term. 

3.1.5 Semiclassical T-term: zone A2 

Tauberian estimate. 

Tauberian estimate for cut-off expression is rather simple: 

5 1 
Cfih-^lo^liV X 'li^ X 7o7i7l«2 - '72~^"2 

which nicely sums to C/i without logarithm due to the same positive eigen- 
value arguments as before; for ^^-cut-off with 9 > iih we get the same albeit 
with 7y defined by the same formula albeit with {wj — IJiJ.hsJ'^) replaced by 
9sj'^ where Sy~^ means the scale; and this should be multiplied by 9/{iih). 
The result nicely sums to Ch~^. This is what was required. 

Magnetic Weyl expression. 

Now we will get the same answer albeit C/x"^ term will be supplemented by 
Cjj,~2h which in cut-off sum adds C/i'^h x ij,h~^ < Ch~^. 
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We can use the standard approach, with an error CiJ,~2h x 9/(fih) x 
IJ,h~'^ X COfi^^h^^ which means that we can take 9 — fi^h which is sufficient 

3 2 

to deal with with 6 > C^^h; in particular, for > hi we are done. But 
for 9 > iih}-~^ we can apply the weak magnetic field approach, which is 
sufficient. So we arrive to inequality 

(3.1.32) I j j 0(x)(e^(x,x,T)-P^,^(M/(x) + T))c/xc/T| < Ch-^ 
and therefore 

Proposition 3.1.8. (i) As B < contribution of zone Xi to expression 

(3.1.33) 1' 1 0(x)(e^(x,x,T)-P^,^(l/l/(x) + T))c/xc/T 

does not exceed Cmax {{Z + B)\z\,ZiBl); 



(a) As B < Z contribution of zone X2 to expression (3.1.33) does not exceed 
CZ\ 



3.1.6 Mollification errors 

We need to estimate 



(3.1.34) j 0(x)(p^(l/l/(x) + r)) - P's{Wr{x) + T))) dx, 

(3.1.35) J 0(x)(Pe(l/l/(x) + T))-PB(l/l/J'^(x) + r)))c/x, 



(3.1.36) D(0(x)(P^(l/l/(x) + r)) - P^(l/1/J'(x) + r))), 

0(x)(P^(l/l/(x) + r)) - P's{Wr{x) + T)))) 

and 

(3.1.37) ||0(x)V(l/l/(x)-l/l/J^(x))|p. 

We start from local versions (so in fact we dealing with w and wJ^). 
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Obviously after all rescalings /i* = h'jQ ^ and therefore e = hi = 

3 ^3 

^^(To^'^Ti ^72"^)"^ where we set 5 = but we will show that we have a reserve 
to set it as 5 > if we want to estimate (3.1.34) by h and (3.1.35)-(3.1.37) 
by h^. 

We claim that 

(3.1.38) \w- wY\ < Cc; := C^TylToTiT^^) 
and 

(3.1.39) |V(^ - wy)\ < C^i CM{^l^hlef. 

Really it follows from equation (3.1.2). 

Then the contribution of Q;2-element to (3.1.35) does not exceed C^ea\ 
as measure of zone of Q;2-element where w ^ wj^ is 0{£a\). One can see 

easily that <,e = 0{h^) and therefore C<^Ea2 — OipAa^ and the summation 
over Q;2-partition of 1-clcmcnt returns Oijii). 

Modulo above calculations the contribution of a2-element to (3.1.34) 

1 3 -i 

does not exceed Cf3q'ieoL\. One can check easily that qe = 0(^1272 ^) and 
1 7 5 

therefore C<j2£a^ = 0(/i3a^7|) and the summation over a2-partition of 
cci-element returns 0{h^aX) and then the summation over cci-partition of 
1-element returns 0(/i2). 

Similarly, expression (3.1.36) with = does not exceed C^e2<^2 — 
Clv'a\ and the summation over Q;2-partition of 1-element returns 0{t\^). 
However we need to consider disjoint pairs of Q;2-elements belonging to given 
ai-element and their contribution does not exceed 

J 72x'72~y \x - dxdy < Ch^a\ 

3 

and then summation over ai-partition of 1-element returns 0{h^). We need 
also to consider disjoint pairs of ai-elements belonging to given 1-element 
and their contribution does not exceed Ctr' / |x — y\~'^ dxdy = 0{h^). 

Finally, contribution of Q;2-element to (3.1.37) does not exceed C<;lea2 and 

one can check easily that this does not exceed CaYil^^ ^^e summation 
over Q;2-partition of cci-element returns Ca\h^] then summation over 1- 
partition of 1-element returns 0[hi). 

So, scaled versions of (3.1.34) and (3.1.35)-(3.1.37) do not exceed Ch 
and Ch?- respectively. Then original versions of (3.1.34), (3.1.35), (3.1.35), 
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and (3.1.37) do not exceed CCH^ x (C^)"^ = CCH^, C^i^ x (C^)-^ = CCX 
CC^f.^ X (C^)-2 = CCP, and CC'£ x {C£)-^ = C^r^ < CCH\ 
Leaving easy details to the reader we arrive to 

Proposition 3.1.9. (i) Contribution of zone X2 to the mollification error 

2 

(3.1.34) does not exceed CZ3; 

(a) Contribution of zone X2 to the mollification error (3.1.35) does not 
exceed Cmax(Z3, B^Zi); 

(Hi) Contributions of zone X2 to the mollification errors s (3. 1.36) and 
(3.1.36) do not exceed CZ^ . 

and 

Proposition 3.1.10. Let B <Z. Then 

(i) Contribution of zone X2 to the mollification error (3.1.34) does not 
exceed CB^Zl'^; 

(a) Contributions of zone X2 to the mollification errors (3.1.36)-(3.1.36) 

do not exceed CB^Zi^^ . 

Remark 3.1.11. Obviously scaled e = 7o7i727(^/7o7i 72)^"*^ — h.^~^ ^ 
(/x~-^/i)2~'^i which makes reduction possible. 

3.1.7 Zone Xj, 

4 

Zone As defined by ixh > Cq, h < 1, {^(x) < €q?} appears only as 6 > 
here Wj^ is smooth and no mollification is necessary. Further here the 
canonical form contains only one number j = and |D°l/l/| < QC^^"'"' and 
1/1/ xC^. 

Therefore we have non-dcgcncracy condition fulfilled and applying the 
standard theory we conclude that in the scaled version contribution of 6(0, 1) 
to the semiclassical errors in N- and T-terms and into D-term are Ciih~^, 
C/i and C^^h~^ respectively. 

In the scaled version they become CBf, CB£C < CBZ2£2 and CBH^ and 

after summation (where for D-term we need to consider mixed contribution 

of different layers) we arrive to the same expressions calculated as ^ = r = 

21 12 43 43 

6~5Z5 i.e. C65Z5, CSsZs and CBsZs respectively. Thus we have proven 
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Proposition 3.1.12. Let < B < Z^. Then 

2 1 

(i) Contribution of zone to the H-error does not exceed CZsBs; 

(ii) Contributions of zone A3 to the T-error and D-term do not exceed 
CZIBi. 



4 Analysis in the boundary strip 
4.1 Properties of as /V = Z 

To finish our analysis we need to get the same estimates as before in the 
strip 

(4.1.1) y -^{W + v< eZ?-\ er < £(x) < c?} 

with ? X B~^ for B < Z^ and r x B^lzl for B > Z^. Let us rescale 
X i-> x' = xr~^. Then 

(a) In the case 6 < we need to rescale then w{x') — B~^Wj^{x'?) and 
take h = 6~4 < 1^ /j, — B^ > 1, nh = 1. One can treat the case 6 < 1 
exactly as in paper [IvrlO]. 

(b) On the other hand, for B > one should set w{x') — rN~^W^^{x'?) 
and h = (/V?)-i = {BN-^)-5 <!,// = fi/V-iri = fii/V-i > 1, i^h = 

bIn-I > 1 (//^x 1 iff e < cnI, nx 1 iff e X n^). 

We will use now only rescaled coordinates unless the opposite is specified. 
Then in our zone 

(4.1.2) Aw = Kw^, K — 12, w ^ 6 — i^C~^ kl ~^ 
with 

' e for 6 < Nh 

N?-^ = bInI for B> nI 



(4.1.3) 



where one can always get k = 12 after rescaling w i-> 144k ^w. 

And right now we consider only case B > (Z — N)i. 
We skip bar over w for simplicity here. 
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Proposition 4.1.1. Let Z = N. Then in y after rescaling 

(4.1.4) \D"w\ < Cai/i/7"l"l Va 

with the scaling function 7 = m/? and 

(4.1.5) |Vm/^| < 1 + Cm/' 
with some constant C and exponent t > 0. 

Proof. Rescaling x 1— )■ x?~^ we get an equations (4.1.2)^ = (4.1.2) with = 0. 
We know that 1/1/ = for £(x) > cr; so after rescaling w — Q for ^{x) > c. 
On the other hand, m/ x 1 as £{x) < e (uniformly with respect to all the 
parameters) . 

Let us consider solution of the equation 

(4.1.6) Aw, = 12w^ 

in Q = {w < e, £ < c} with boundary conditions = w at dO.; s > |. 

Note first that > 0. Really, is the solution of the variational 
problem to minimize ||Vm/|P + 24(s + w^^dx and one makes this 

functional only less replacing 1/1/ by 

Further, the standard maximum principle arguments show that \ as 
s\'^^\ Obviously 1/14 \ 1/1/ and i/i/j — > 1/1/ in in {x : w{x) > 0} as s \ |. 

We claim that 

(4.1.7) Ws e ^^^+2 

To prove (4.1.7) note first that w G uniformly with respect to all 

the parameters for any 6' > 0. Then G and then (4.1.6) yields that 

1 

Ws G "^^^^"^ as soon as s — 5 ^ Z. Then since w, > we get \Vws\ < cwi 

1 _ 3 

and so 1/1// G ^^~2. Then equation (4.1.6) again yields that w, G Now 
we need more subtle arguments. 
First, for |y| = 1 

(4.1.8) ws{x + ty) = ws{x) + t{Vws), ■ y + \{V^Ws),{y)t^ + 0{e). 

If Aw/ = fi{wi) in Q, fi{w) sls w and /i(M/) > /i(w) then A(wi — W2) > as 
wi > 1^2 and then w\ — does not reach maximum inside Q. 
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Then the lowest eigenvalue of V^m/s at x should be greater than —Cwi. 
Really, otherwise we can take y as the corresponding eigenvector and t 
with |t| = and with a sign making second term non-positive and get 
Ws{x + ty) < 0. 

This lower estimate for eigenvalues of V^m/s and equation (4.1.6) yield 

1 2 

that |V^M/s| < Cws ■ But then |Vm/s| < Cwi. Really, otherwise picking 
y = |Vm/s|~'^Vm/s with |t| = e|VM/sh appropriate sign we would get 

Ws{x + ty) < 0. 

1 

These estimates yield that Ws{x') x Ws{x) in 6(x, 7(x)) with 7(x) = ewi . 
Then e 'rfh In fact, let us consider f = w^~^\/w and \ f{x) — f{x')\. Let 
us consider first |x — x'\ > |(7(x) + 7(x')); since \f{x)\ < 7(x)2 at each 
point we get that \f{x) — f{x') < \x — x'\^. 

On the other hand, for |x — x'| < |(7(x) + we conclude that 

7(x) X 7(x') and \f{x) — f{x')\ < |Vf | ■ |x — x'\ < |x — x'|^ because 
\Vf \ < |VVs|m/5*-^ + |Vm/sP<~^ < C7^-\ 

Therefore e and equation (4.1.6) yields that i/i/j G 

In the next round we assume that w e '^^+^~'' with some S e (0, 1). 
Then 

(4.1.9) Ws{x+ty)< 

w,{x) + t(Vi/i/s)x • y + J(V2i/i/s)x(y)t' + ^(W,),(y)t^ + C\t\P 

Z 

with p = min(4, 4 + s — 6). 

We claim now that the lowest possible eigenvalue q of (V^i/i/j)^ is greater 

than —Cws''~^^^'^. Really, otherwise let us pick up y as the corresponding 
eigenvector, t with |t| = e|^^|-'^/(''~^^ and with a sign making expression 

t{Vws),-y+^t\V^WsUy) 

non-positive and get Ws{x + ty) < again. Now equation (4.1.6) yields that 
inequality 

(4.1.10) ^ I W,| < Ci/i/i"-'^)/" 
holds with k = 2. 
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Further, we claim that this inequahty holds with k = 1,3. Indeed, if one 
or both of these inequalities are violated then let us take corresponding y 
and t with 

\t\ = ellVwsl'^^"-'^ + \V'w,{y)\'/^P-'^^ 
(calculated on y); replacing e by 2e if necessary we get 

|t(Vi/i/s)x-y + ^t^(Ws)x(y)| > eo|t(Vi/i/s)x-y| + |^t^(Ws)x(y)| 

and choosing an appropriate sign of t we get w{x + ty) < 0. 

Therefore inequalities (4.1.10)j_3 hold. The same arguments as above 

with 7 = Ws^'^ lead us to G ^''^ and then equation (4.1.6) yields that 
i/i/s € 't^'^^^^. So, now we came back with 6 replaced by 5' = 2 + s — ps and 
one can see easily that if 5 > s then 6' — s + {2 — 4s) + {6 — s)s and after 
few repeats S < s. Then we get (4.1.7). Unfortunately, constants depend on 
s due to the fact that Am/ e fails to yield w e^^. 

Now we are going to finish the proof of (4.1.4). Let us consider Ws again 
and let 7 = 75 5 = Ws^^^~^\ Due to the previous inequalities 7 G We 
claim that IV7I is bounded uniformly with respect to s, S. Note first that 
^^4-5 ^ ^(4-5)s iinplies that 

(4.1.11) a|V7p + 67A7 = 7^^ 

with a = ^(4 - S){3 - S), b ^ ^(4 - S), and a = 4s - 2 + (1 - s)5. 
Let ^0 = IV7P; obviously t/j is uniformly bounded at dO.. Let us consider 
maximum of xjj reached inside Q. At the point of maximum 

(4.1.12) ^7,,.,.7,,. = 
and 

^ = y: 7ix. + E 7. (At), = E ^l, + E ^>^-- (7-^ - ai V7n) 

iJ i ij i 

due to (4.1.11) and due to (4.1.12) this expression is equal to 

E^x,x, - b-'r'\V^\'{r - a|V7p) + 6-V7^-^|V7p 

iJ 
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and therefore at an inner point of minimum a|V7p < 7*^. So, IV7I < C is 
proven and for s \ |, 5 \ we get that | Vi/v^ | < C. 

Let us pick 7(x) = e'w^(x); then |Viv| < | and w[x) x w[x) in 
6(x, 7(x)). This and equation (4.1.6) easily yield (4.1.4). 

To prove inequality (4.1.5) let us consider Ws again and let us take 
now ip = IV7P — M^'^^ with t > 0; obviously ip is non-positive at dO. for 
sufficiently large M. Let us consider maximum of ip reached inside Q. At 
the point of maximum 

(4.L12)' Yl ^/x.-.^-'- - = 

and the same arguments as before (plus inequality IV7I < Q) show that 
at an inner point of maximum a|V7p < 7^^ + CtM^^^ where C does not 
depend on M and small t > 0. Then at this point ijj <1 for small enough 
t > and as s ^ I and 5 ^ we get (4.1.5). □ 

The following statement heavily uses estimate (4.1.5): 
Proposition 4.1.2. The following estimate holds 

(4.1.13) D(7-^+^ 7-^+^) < 

with some constant C which does not depend on s e (0, 1) where we set 
:= wl^"^^'^ (i.e. itisO as w <0). 

Proof. As in the notations of the proof of proposition 4.1.1 5 = and s = | 
we have (4.1.11) with a = 1, 6 = 3 and a — 0: 

(4.1.14) ^7A7+ |V7|^ = 1. 
Then 

7-'+^ = 7-'+lV7|^ + ^7^A7 = (1 - |)7-'+lV7|^ + ^(^^^7'+^ 
and 

D(7-'+^ 7-'+^) < (1 - |)D(7-^+lV7|^7-'^^) + C < 

(1 - |)D(7-^+^7-^+^) + CD(7-^+^+^ 7-^+^) + C 
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due to (4.1.5) and this yields 

D(7~''''-7~'+') < Cs-'D(7-l+^+^7-l+^+^) + Cs-\ 
Substituting s + mt instead of s, < m < Ct~^ we recover (4.1.13). □ 

4.2 Semiclassical analysis in the strip y for 
N>Z. 

It is really easy to construct the proper potential in this case: we just take 

(4.2.1) w, = w(t)„ (P, = fiwe-"") 

with f e ^o°°((|, oo)), < f <1, f{t) = 1 for t > 1. Note that due to 
(4.1.2) 

D(7-V.,7-V.) < Ce-''D{Y-\Y-') < Cs-h-'\ 
D(l - 0„ 1 - 0,) < Ce'-''D{Y-\ Y-') < Cs-'e'-'^; 

then minimizing with respect to s (= | log£|~'^) the right-hand expression 
we conclude that 

(4.2.2) D(7-V., 7"V.) + £~'D(1 -<f>,,l- <f>,) < C(l + | log£|)' 
and therefore 

(4.2.3) J j-^cf>,dx + e-^ J{l-<l),)dx<C{l + \\ogs\). 

Remark 4^.2. 1. (i) Recall that all these integrals are taken over domain 
{x : w[x) > 0}. To avoid possible troubles we pick e = and set in the 

4 

zone {x : w{x) < Coh^} 



(4.2.1)' 



7(x) = dist(x,{M/>2Con3}), 

—■~f^4>'g for 7 < e, 
—e^ for 7 > £ 



with (j)'^ = {{-js ^) and then in the complemental domain {x : w{x) < —<;^} 
with <, — and <J7 = 7^ > h; 
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(ii) Further, for e — h3 ^ with sufficiently small exponent 5 > it does not 
break estimate for mollification error in T-term; 

(iii) Furthermore, for t > £ 

mes({x : 7(x) < t}) < Ct^e'^ mes({x : 7(x) < ee}) 

and therefore 

h"" j 7-1-"'^-"' dx < Ce-^ mes({x : 7(x) < te]) < CL := C(l + | log h\) 

for sufficiently large m. 

Using these estimates and the last remark we can prove easily 
Proposition 4.2.2. Let N > Z. Then 

(i) Contribution of y U with X4. := {x : w{x) = 0} to mollification 
and approximation errors in N-term do not exceed CToe^{l + \ log£|) and 
/?o(l + I log^l) respectively with 

(4.2.4)1 Tq^b'^, Rq^bK t^bI, r = b-^ 

as B <Z'3 

and 

(4.2.4) 2 Tq^z, /?o = fi^zi, t^zIb-5, r = zIb*^ 

as Zl <B < 

(ii) Contribution ofyuX/i to mollification and semiclassical D-terms do 
not exceed CTe^{l + | log£|) and R{1 + \ log£|)^ respectively; 

(iii) Contribution o/3^U to both mollification and approximation errors 
in T-term do not exceed CTe^{l + \ log£|) and CR respectively. 

Proof. Really, estimates for mollification errors and terms immediately follow 
from the inequality 

(4.2.5) mes({x : w[x) < £^}) < Ce{l + | log£|) 

which is due to (4.2.3). 

Let us consider semiclassical errors and terms. 
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(i) Let us consider A/-term first. The contribution of eacli ball 6(x,7(x)) 
to the semiclassical error docs not exceed Cfih~^'~f^ x CB?^ and the to- 
tal contribution does not exceed CRq J dx < C/?o(l + | log£|) where 
Ro = Br^; 

(ii) Consider semiclassical D-term. The similar arguments with the analysis 
of disjoint balls of different types and with analysis of the intersecting balls 
(of the same type) lead us to the proper estimate of the contribution of 3^4UA'4 
to semiclassical D-term: namely, it does not exceed CRq?~^(1 + | log£|)^ 
where Rq?^^ x /?; 

(iii) Consider T-term. Let us consider all possible balls and their contribu- 
tions. The contribution of each ball 6(x,7(x)) to the semiclassical error 
does not exceed C^^fiq^^ x CB(^?qj^ and the total contribution does not 
exceed CR J q-f'^ dx x CR where R = 8(^7. □ 

Then picking appropriate e = we conclude that 
Corollary 4.2.3. Let N>Z. Then 

(i) Contributions of y U X/^ to all errors in H-terms do not exceed CRqL 
with {l + \\og BZ-^\); 

(ii) Contribution of y U X4 to all D-terms do not exceed CRL^; 

(iii) Contribution o/3^ U A4 to all errors in T -terms do not exceed CRL. 

Remark 4-2. 4- Could we get rid off logarithmic factors i.e. make /. = 1 as 
in the case M = 1? 

(i) With the mollification errors we need to replace (4.2.5) by 
(4.2.6) mes({x : w{x) < e"^}) < Ce; 

(ii) With the semiclassical terms our arguments here are insufficient even if 
we established (4.2.6); we need extra propagation arguments in the direction 
of decaying w along magnetic lines-exactly as in the case M = 1. Surely 
there could be points where such arguments do not work; f.e. consider 
M = 2 and nuclei so that |yi — 72] is slightly less than ?i + ?2 where ri 2 are 
precise radii of support. Then w reaches its minimum at y. 

So, we need to prove that the measure of such points is sufficiently small 
(f.e. less than C| \ogBZ-^\-^). 
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Unfortunately I do not know how to make the above remark work and I 
suggest 

Problem 4.2.5. Follow through the discussed plan. For M = 2 it could be 
easier due to rotational symmetry. 

4.3 Analysis in the boundary strip y. 

Construction of the potential for N < Z 

Now let us consider N < Z. 

Remark 4-3.1. (i) We can replace < by = in the definition of Wj^ 
without increasing any of errors as long as Br^u2[i//(^)i < Rq; 

(ii) As 6 < Zt this translates into |z/|t < zt i.e. (Z — A/)+ < Z^B~^; 

(iii) On the other hand, as Zi < < Z'^ it translates into < B^/^ i.e. 
{Z- N)+ < Z'sBrs. 

Therefore let us consider a case not covered by remark 4.3.1. Let 
9 = < 0. 

Proposition 4.3.2. Consider dependence of W^^ = Wj^^-^{x) on v. Then 

(i) V\/B(v)i^) + ^ non- decreasing with respect to v at each point x; 

(ii) y^B(u)i^) '^^ non-increasing with respect to v at each point x; 

(iii) In particular Wj^^^{x) + v l/l/J(o)(x) and Wl^^^{x) \ l/l/J(o)(^) «^ 
each point x as u ^ —0. 

Proof, (i) Consider Wj = ^b(v) + ^7' > ui > U2. One can prove easily 
that Wj^ — V is a. continuous function and since l/l/i — I/I/2 = i^i — ^2 > as 
i[x) > C we conclude that l/l/i > I/I/2 at each point x (which is exactly our 
statement (i)) unless l/l/i — I/I/2 achieves a negative minimum at some point 

X*. 

(a) Let X* ^ y„; then A(l/l/i - l/l/2)(x*) = P^(lVi) - ^6(1^2) < because 
H/i < W2 at X* and x* cannot be such point. 
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(b) Let X* = ym- Prom Thomas- Fermi equations for l/l/i 2 one can prove 
easily that 



(l/l/i - l/l/2)(x) = (l/l/i - l/l/2)(y„) + U{x - y„)+ 

«:;„|x-y^| 5(1/1/1 - l/l/2)(ym) + 0(|x-y„|2) 

near y^ where Lm{x) is a hnear function and > and therefore if 
(1/1/1 — l/l/2)(ym) < 0, ym cannot be a minimum point either. 

(ii) So, 1/1/1 > 1/1/2 cind therefore l/l/i — I/I/2 is a subharmonic function and 
since l/l/i — I/I/2 = z/i — z/2 as £(x) > C we conclude that l/l/i — I/I/2 < i^i — ^2 

'^e.Ci/i) - ^BXu2) ^^"^^ point. 

(iii) Statement (iii) follows from (i), (ii). □ 

Therefore in the zone {x : 1/1/J(^) > (1 + e)|i/| we can apply the same 

(7, (;) scaling with = 7^ defined for = 0. Really, we know that there 

KU + ^ ^ Klo) ^ and , = 7^. 

Then Thomas-Fermi equation (1.1.3) implies that 

|V"14/Jf,)| < C^V'"' Va 

and then we arrive to the statement (i) in Proposition 4.3.3 below. 

On the other hand, in the zone {x : 1/1/J(^) < (1 — we can apply 

the same arguments but this zone is classically forbidden and we arrive to 

to the statement (ii) below. In both cases ^7 > (where in the latter case 7 
is the distance from x to H/J^f"^) (scaled) and <; — |z/|2 because in both cases 

^" > I'^I^C""^ and 7 > u^^^i so <^7 > I^^I^C^^. 

Then we arrive to the statement (ii) in Proposition 4.3.3 below. 

Proposition 4.3.3. Let either B < and > Zl or < B < Z^ 

I 1 3 1 
and > 62. Then 

(i) Contributions of zone {x : l/l/e'^(x) > (1 -|- eo)]!'!} to semiclassical errors 
in N- and T -terms and into semiclassical D-term do not exceed CRqL, CR 
and CRL? respectively; 

(ii) Contributions of zone {x : W^^^x) < (1 — eo)|i^|} to semiclassical errors 
inH- and T -terms and into semiclassical D-term do not exceed CRqL, CR 
and CRL^ respectively. 
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Remark 4 -3. 4- (i) Here actually /. = 1 + | log ^| with 

(4.3.1) \u\C' > 

(ii) To have a proper estimate for D-term with i/ replaced by we need 

(4.3.2) (Z - NU < 



{Z-N)+B--^ as {Z - N)l< B < Z-\ 
(Z - /V)+Z-^ as Z5 < S < Z^ 



-, 10 „ 25 „ _ 4 

ZTg-is as e < Z3, 



Z^B^ as Zl < B < Z^ 



,5 



with the right-hand expression larger than Zs. 

Therefore we need to explore zone 

y* := {x : (1 - eo)|i^| < W^'{x) < (1 + eo)|i^|} 

in the framework of Proposition 4.3.3. 

We assume that /i* < 1 where 9 = |j^|C~^ and 

(Z - A/)+^ei5 as e < Z5 



(4.3.3) K = hO'-^ 



{Z-N)^'Z^^bI asZl<B<Z 
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This assumption is a bit more relaxed that the negation of the framework 
of Proposition 4.3.2. 

Proposition 4.3.5. Assume that condition (4.3.2) is violated. Let us con- 
sider B{x, T)) (after reseating x i->- x?"^) with r] := 9^ and scale it x ^ xr]~^. 
Let w := (l/l/J^ + z/)|z/|-^ 

Then there exist scaling functions cto, ... , 0:3 with aj = ^yjaj^i, ccq = Jo0~i 
such that 

(4.3.4) |V"M/| < C^o'To"'"' ^ith ^0 = To. 

(4.3.5) 70 X nnax(|M/|?, |Viv|^ |V^m/|^, |V^m/|), 

and after rescaling x x'Jo'^, w wq — w^q'^ 

(4.3.4) , iV^iVol < QqV'"' ^ith ^1 = 7I, 

(4.3. 5) ;^ 71 X maxdwol^, |Vi/i/o|5, IV^wqI), 
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and after rescaling x i->- x^^ ^, wq wi = wq^^ ^ 

(4.3.4) 2 |V"l/l/l| < C«^2'72"'"' with ^2 = 72- 

(4.3.5) 2 72 X max(|M/|^, |Vm/|), 
and after rescaling x ^ X'j2^, 1/1/1 i—)- 1/1/2 = 1/1/1^2^^ 

(4.3.4) 3 |V'^iV2| < C^|73~'"' with ^3 = 73~'- 

(4.3.5) 3 73 - |>^|- 

Proof. Proof easily follows the arguments of the beginning of the Section 3; 
easy details we leave to the reader. □ 

If needed we replace 7/ by max(7j, 7y) with 

21 2 
7i — CHIq, 72 := Chl -^, 73 := C^^^2. 

_2 _1 

^*,o ^*7o^ . ■= ^*,o7i ^ ^*,2 := ^*,i72 ^■ 

We leave to the reader to prove by our standard methods 

Proposition 4.3.6. In the framework of Proposition 4-3.5 consider a,- 
element 0/. Then one can apply e; -mollification to potential 1/1// with e-, — 7; 
as 7/ < 27/ such that 

(i) Contributions of (pj-element to both semAclassical and approximation 
errors in U-term do not exceed CB?^a\, CBT^a\, CB?^a\{l + | \ogh^^i\) and 
CBr^al{l + I log^*,i|)^ for i = 3, 2, 1, respectively; 

(a) Contribution of {(p;, (j);) -pair to both semiclassical and approximation 
in D-terms do not exceed CB'^T^al, CB^T^a^, CB^T^a^l + | log/i*,i|)^, and 
CB^r^aKl + I log^*,i|)'* for i = 3, 2, 1, respectively; 

(Hi) Contributions of (pi-element to both semiclassical and approximation 
errors in J-term do not exceed CB^?^a\, CB^r^al, CB'^?^a\, and CB^T'^Oq 
for / = 3, 2, 1, respectively. 

Then we arrive to 
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Proposition 4.3.7. In the framework of Proposition 4-3.5 there exists po- 
tential such that 

(i) Contributions ofy* to both semiclassical and approximation errors for 
U-term do not exceed CBr^J with 

(4.3.6) J:= f ao{x)-\l + \\ogK%^\fdx 

where ?~^y* is y* scaled by x i-^ xr~^) and ao — 7o^~^; 

(a) Contributions ofy* x y* to both semiclassical and approximation D- 
terms do not exceed CB^T^J^; 

(Hi) Contributions ofy* to both semiclassical and approximation errors for 
T-term do not exceed CB^r^ . 

Remark 4-3.8. Obviously 

(4.3.7) J < a max [ 7o(x)-'(l + | logn,7o-^|)' dx 

JB(0,1) 

where this integral is taken over B{x, CqO^) rescaled which in turn does not 
_ 1 

exceed Ch* ^ which translates into 

(4.3.8) (Z - N)Il X 
for errors in N-term and 

(4.3.9) (Z - N)IL^ X 

for D-terms. 



as B <Zi , 

Z^oB^ as Z-3 < B <Z^ 

B2i as B < Z\ 

Z'^B^ as Z'i < B <Z^ 
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4.4 Conclusion 

Let us recall that 



(4.4.1) 



(4.4.2) 



and therefore 




as B<{Z- N)l, 

as {Z-N)l<B< Zt, 

as < e < Z^ 

as B< {Z-N)l, 

as B <Z^, 

as Z'^ < B < Z^ 



(4.4.3) K = h9—^ = < 



and 



(4.4.4) 



{z- 


N)7 


as 


e < (z - N)i, 


' (z- 


A/)7e^ 


as 


{Z-N)l<B< Zt 


Az- 


A/);^z^e^ 


as 


Zl <B<Z^ 


[z- 


A/)7e 


as 


B < (Z- A/)|, 


\t 


/v)|e5^ 


as 


(z-/v)| < e < Z5 




/\/)|z^e^ 


as 


zt < e < z^ 



Therefore we arrive 
Proposition 4.4.1. Let M >2. Then for a constructed above potential W 
(i) Total semiclassical and approximation errors in H-term do not exceed 



Z-3 +{Z - N)_^'B 



(4.4.5) C < 



Z3 + B2L + {Z - N)IB^L* 
zlBh + {Z - N)lzi^B^L* 

where here and below /. = (1 + | log BZ^\), /.* = (1 + | log 6*1); 



as B<{Z- N)l, 
as {Z-N)l<B< Zl, 
as Z'i < B < Z^ 
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(ii) Both semiclassical and approximation D-terms do not exceed 



Z'^ + B^{Z - N)+' asB<{Z-N)l, 



(4.4.6) C < 



Z-3 + B-^L + {Z - N)lB2iL* as {Z - N)l< B < Z-3, 
zIbIl + {Z - N)lz'2Bh* asZl<B<Z^; 

(Hi) Total approximation error in T-term does not exceed 



(4.4.7) C 



Zi + e?/. as B < Z5, 



zIbIl as Z-3 < B < Z^; 

(iv) Total semiclassical error in T-term does not exceed 
(4.4.8) (4.4.7) + CZifi^ + CZia"^ 

provided a > Z~^; as a < Z"^ the last term should be replaced by CZ^. 
Also we arrive to 

Proposition 4.4.2. Let M > 2, B < Z and a > Z~\ . Then for a con- 
structed above potential W 

(i) Total semiclassical and approximation errors in N-term do not exceed 
CZl{{BZ-y + (aZ5)-5 ^ z-^^^. 

(ii) Both semiclassical and approximation D-terms do not exceed and total 
semiclassical and approximation errors in T-term do not exceed CZ^ (^[BZ~^Y-\- 
{aZl)-' + Z-'). 



5 Ground state energy 
5.1 Lower estimates 

Now the lower estimates are already proven: in virtue of the analysis given 
in the previous paper [IvrlO] these estimates are due to estimates of the 
semiclassical and approximation errors for X — v. 
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Proposition 5.1.1. Let M = 1. Then 
(i) The following estimate holds: 

(5.1.1) E"^'= > + {lr{HA,w + J Pb( 1/1/"^'= + u) dx^ - CR; 

(a) The following estimate holds: 

(5.1.2) E'^^ > 8'^^ + Scott -CR-CZIb'^; 

(Hi) As B <Z 

(5.1.3) E'^^ > E^^ + Scott + Dirac + Schwinger - CZi (Z"* + [BZ-^)^) . 

Proposition 5.1.2. Let M > 2. Then 
(i) The following estimate holds: 

(5.1.4) E"^"^ > S"^ + {lr{HA^w - + J PsiW^'^ + u) dx) - (4.4.7); 
(a) The following estimate holds: 

(5.1.5) E'^^ > S^^ + Scott - (4.4.7) - CZtfii - CZia'^ 

provided a > Z~^ and B < Z^; otherwise we can skip Scott and replace the 
last term by CZ^; 

(Hi) As B < Z and a > Z~\ 

(5.1.6) E'^^ > £'^^ + Scott + Dirac + Schwinger- 

CZl{Z-^ + {BZ-y + {aZ's)-^). 
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5.2 Upper estimate: general scheme 

Recall that for an upper estimate in addition to the trace we need to estimate 
also \Xi\i — u] and 

(5.2.1) |Aa/ - z/| ■ |N(H^,M/) - A/| 
and also three D-terms: two of them semiclassical 

(5.2.2) i 2 D(e(x, x, A) - P's{W{x) + A), e(x, x, A) - P'b{W{x) + A)) 

with \~v and A = A/v and also 
(5.2.3) 

D (P^( H/(x) + A/v) - P'b{^{A ^v),P's{ H/(x) + Aw) - l/l/(x) + i/)) ; 
our tool will be semiclassical estimates for two semiclassical N-terms 

(5.2.4) i2 y (e(x,x,A)-P^(l/l/(x) + A))cyx 

with \ = v and A = A/v and also estimate from below for the third N-term 

(5.2.5) I / + A/v) - P'b{ + y)) dx. 

5.3 Upper estimate as M = 1 
5.3.1 Estimate for |A/v — 

As in paper [IvrlO] we have two cases: in the first case u is small enough so 
we construct W'^^ with u = and estimate |A/\/| and in the second case u 
we prove that Aa/ x i/ and estimate |Aa/ — 

Proposition 5.3.1. Let M = 1, B < Z^. 

(i) Let 

(5.3.1) {Z-N)+ < K:= Comax(zi, zifi^). 
and let us construct W as if u = i.e. N = Z. Then 

(5.3.2) |A/v| < Cimax(zi, ei). 
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(ii) Let 

(5.3.3) (Z-/V)+ > K = Comax(zi, Zie^) 

with sufficiently large Cq. Then A/v x and 

(5.3.4) |Aw-i/| < Cimax(zi, fi^)|i/|?. 

Proof, (i) One can see easily that 

4 

(5.3.5) In the framework of (i) as 6 > (Z — N)^ expression (5.2.5) is 

X B|A/v|5 X (^-^y~r'^ ^ |Aa/|^ min(l, S~razi) 



where (IAa/I/S^)^ is a width of the zone where W < —Aw and the se- 
lected factor is the volume of this zone; this expression should be less than 
Cnnax(Z3, Z565) which is exactly an error estimate in the semiclassical 
expression for N. Thus 

(5.3.6) |Aa/|^ min(l, S'^zi) < Cmax(zi, zifi^) 

4 4 

where everywhere the first and the second cases are as 6 < Zs and Zs < 
B < respectively. The last inequahty is equivalent to (25-6-13). 

— 3 2 

As 6 < (Z — /V)^ inequality (5.3.6) is replaced by \Xn\^ < CZs which 

4 

coincides with (5.3.6) with B reset to (Z — N)\. and also with the same 
inequality derived for 6 = in the previous Chapter; therefore (5.3.2) holds 
in this case as well. 

(ii) One can prove easily that 

(5.3.7) if (5.3.3) is fulfilled, and expression (5.2.5) does not exceed a semi- 
classical error then A/v x t/ and, furthermore, expression (5.2.5) is 

(5.3.8) xe|A,v-i^| J Pb{W + ly) dx >^ B\Xn - iy\ J {W + ly)^^ dx 
which for B > {Z - N)l is 

(5.3.9) X B?^\Xn - u\ ■ (-^) ' X |A/v - i/| • min(l, fi'^zi) 



76 



and this should be less than Cmax(Z3, Zsfis) which implies (5.3.4). 

4 

On the other hand, as 6 < (Z — A/)^ the right-hand expression of (5.3.8) is 

X |A/v — z/|r X |A/v — v\{Z — N)_^_^ and this should be less than CZs which 
implies (5.3.4) in this case as well. □ 

Proposition 5.3.1 immediately implies 
Corollciry 5.3.2. In the framework of proposition (i), (ii), 

(5.3.10) \Xn-i^\ ■ N([AA/,t/]) < Cmax(Zi,zist) 

where N(A/v, i') is the number of (non-zero) eigenvalues on interval [A/v, 
or \v, A/v]. 

5.3.2 Estimate for D-terms 

Proposition 5.3.3. In the framework of Proposition (i), (ii) 

(5.3.11) D(e(x,x,A)-P^(l/l/ + A), e(x, x, A) - P^(l/1/ + A)) 

with X = 1/, Xi\i respectively and 

(5.3.12) D{P'b{W + u)- P'b{W + X), P'b{W + p) - P'^iW + X„)) 

5 3 4 

with X — Xi\jdo not exceed Cmax(Z3, Z^Bi). 

Proof. Term (5.3.11) with X = v has been estimated this way and the same 
estimate for this term with A = A/v is proven in the same way; we leave easy 
details to the reader; 

Term (5.3.12) is estimated using proposition 5.3.1; again we leave easy 
details to the reader. □ 

Remark 5.3.4- Let B < Z . Then in (5.3.1)-(5.3.4) and therefore also in 
(5.3.8) and in proposition 5.3.3 one can replace Co and C by CqS and Ce 
respectively with the small parameter e: max(Z~'', {81"^)^^ < £ < 1- 
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5.3.3 Summary 

Then following the scheme of subsection 24.4.4 of [Ivrll] we arrive to upper 
estimates in Theorem 5.3.5 below (lower estimates have been proven in 
Proposition 5.1.1); we also arrive to Theorem 5.3.6 below: 

Theorem 5.3.5. Let M = 1, B < Z^. Then 

(i) The following estimate holds: 

(5.3.13) E"^'' < £'^^ + {lr{HA,w - + j Ps( l/l/^''(x) + u) dx^ + CR 

/ 5 3 4 \ 

with R = max(Z3, Z^B^y, 

(a) The following estimate holds: 

(5.3.14) E'^^ < S^^ + Scott + CR - CZl bI, 

(ill) As B < Z 

(5.3.15) E^^ > £'^^ + Scott + Dirac + Schwinger + CZi (Z^^ + {BZ-^) . 

Theorem 5.3.6. Let M = 1, B < Z^ 

(i) The following estimate holds: 

(5.3.16) D(p^-p7,pw 

(ii) As B < Z 

(5.3.17) D{p^-pY,p^-pY)<CZl{Z-' + {BZ-y). 

5.4 Upper estimate as M > 2 
5.4.1 Estimate for |A/v — i^l 

Again we need to consider two cases: almost neutral molecules (systems) 
when (Z — A/)+ < CqK with K slightly redefined below and we can set u = 
in the definition of Thomas- Fermi potential and establish estimate for |A/v| 
(and for optimal u we have the same estimate for both and A/v) and not 
almost neutral molecules (systems) when (Z — A/)+ > CqK and we can prove 
that |Aa/| X \i>\ and estimate |Aa/ — 



Then 

- p7) < CR- 
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Proposition 5.4.1. Let M >2, B < and condition (0.4.1) be fulfilled, 
(i) Let 



,2 



73 + 62/. as B <Zi, 

Z'^B^^L as Z'^ <B <Z^ 



(5.4.1) (Z -N)+<K := Q 
and let us construct W as if u = i.e. N = Z. Then 

(5.4.2) \Xn\ < Ci 



Z9 + Bn3 as B < Z3, 

Z^^B^^L^ as Z-i < B < Z^ 



recall that /. = | log BZ '^\; 
(ii) Let 

(5.4.3) (Z -N)+>K 

with sufficiently large Q. Then Xi\i >i 1/ and moreover 

(5.4.4) |Aa/ - i/| < Cmax(zi, ei/.i)|i/|^; 
recall that 

{1 as S < (Z - /V)!, 

|iog((z- /v)+/e^| as (z-/v)| < e < zt, 
|log(Z- /V)+/6tzi|) asZt<6<Z^ 

Proof. We will apply arguments slightly more sophisticated than the obvious 
ones. These better arguments will allow us to derive slightly better estimates 
for \XN-jy\ as {Z-N)+ > CK. 

Recall that estimates for \Xf\j~u\ are derived by comparison of expression 

(5.2.5) and scmiclassical errors for the number of eigenvalues below X = v and 
A = Aft/: expression (5.2.5) should be less than the sum of these semiclassical 
errors. 

Consider contribution of each ball 

(5.4.6) B{x, ^{x)) c 3^ = {x : min |x - y„| > er} 

to semiclassical errors as A = and A = A/v and compare it with its 
contribution to (5.2.5). 
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(a) Each ball contributes no more than CBe^ to the first error (with X = u) 
where due to our choice ^ 1- 

(b) Further, each ball with ( > Ci\Xi\i — contributes no more than CBi^, 
but each ball with ^ < Ci\\i\i—u\2 contributes no more than CBi'^~^\Xi\j—iy\^/^ 
to the second error (with A = A/v); here 5 = | is due to rescaling. 

(c) Meanwhile each ball with ( > Ci\Xi\i — contributes no less than 
eoB\Xi\j — u\(~^£^, and each ball with ( > Ci|Aa/ — u\2 contributes no less 
than eo|A/v — z^|2£^ to expression (5.2.5) and it is larger than its contributions 
to each semiclassical error (multiplied by C) as long as 

(5.4.7) i 2 ^ I^A/ -i^\> C2Cr\ |A/v -u\> Csr^ 
Obviously in (i), (ii) we can assume that 

(5.4.8) (5.4.2) and (5.4.4) respectively (with C replaced by arbitrarily large 
Cs) are violated. 

(i) Assume first that B > {Z - /\/)|. 

Then in the framework of (i) C = B^i^ with minimal i = B~^ and there- 
fore (5.4.7)i2 are fulfilled for i < C2B^^\Xi\j\. Therefore we need to ac- 
count for semiclassical errors contributed by an inner shell (not exceeding 
C max(Z3, B2)) and by zone y f] {i > C2B^^\Xi\i\}; there ( > Ci|Aa/|2 and 
therefore its contribution docs not exceed CB J ^(x)"-^ dx with integral over 
this zone and it does not exceed CBr^L. 

So, the truncated semiclassical errors do not exceed C max(Zi, Br^L). Mean- 
while expression (5.2.5) is no less than C62r^|A/v|4. Therefore comparing 

these two expressions as 6 < Zs and as Zs < 6 < Z'^ we arrive to (5.4.2). 

4 

Consider remaining case B < (Z — A/)^. Semiclassical arguments remain 
valid while estimate of (5.2.5) from below by eo|A/v|? also could be proven 
easily. 

4 

(ii) Again, assume first that (Z — A/)^ ^ B < Z^. Again, in the calculation of 
the truncated semiclassical errors we integrate over i > C26~^|A/v — where 
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C ^ Q|A/v|2 and therefore its contribution does not exceed CB J £(x) ^ dx 
with integral over this zone and it does not exceed CBr^L\ 

Again, expression (5.3.8) is larger than the expressions afterwards and 
comparing with the semiclassical error estimate we arrive to (5.4.4). 

4 

Consider remaining case B < [Z — N)^. Semiclassical arguments remain 
valid while estimate of (5.2.5) from below by eo|A/v — also could be 
proven easily. □ 

Corollary 5.4.2. In the framework of proposition 5.4-1 |A/v — I'l ■ N([Aa/, i']) 

does not exceed expression (4.4.6). 

5.4.2 Estimate for D-terms under assumption (5.4.1) 

We need to estimate semiclassical error D-tcrm (5.3.11) with A = Aa/ because 
for \ = u we already estimated it, and also we need to estimate a term 
(5.3.12). We start from the latter one. Recall that under assumption (5.4.1) 
we take u = 0. The trivial estimate is based on 

(5.4.9) \P'b{W) - P'b{W + X)\ < C1/1/5|A| + Cfil/l/-^ |A|i, 
leading to 

(5.4.10) J < CD{WK l/l/i)|A|2 + CS2|A|iD(l/l/-^^, 1/1/-^^) 

where here and below J is expression (5.3.12), 6 is a characteristic function of 
domainjx : 7(x) > hi} and we can ignore contribution of {x : 7(x) < hs}. 
Really, contribution of this zone does not exceed a semiclassical error estimate 
CRL^. 

Note that even without assumption (5.4.1) 

(5.4.11) D(l/l/5, 1/1/5) X (e-^; e-izi) as B < z'^ , zl < B < 

respectively and (5.4.2) implies that the first term is much less than R. 
Meanwhile under assumption (5.4.1) 

(5.4.12) D(i/i/-^^, w-'^e) X B-^D{r^e, r^e) x B-^?^D{'y-^e, ^-^e) 

In (i) this would not lead to any improvement. 
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where in the right-hand expression D and 7, 9 are in the scale x 1— )■ xr 
and then 0(7"-^^, 7"-^^) x so the second term in (5.4.10) does not exceed 
CB?^\\i^\il? which due to (5.4.2) does not exceed 

(5.4.13) Cmax(zi, Ziei/.^). 

Consider now term (5.3.11) with A = Aa/. Let us consider zones Q.i '.= 
{x : |A - z/| < Ci-'} and Q2 := {x: \X-iy\> C^i}. 

Remark 5.4-3. Contribution to the term in question of each couple of balls 
contained in Qi x Qi does not exceed estimate for the same term with A = z/; 
really, after rescaling x 1— )■ x/£ and r h- )■ r/C^ we conclude that the difference 
between energy levels does not exceed local semiclassical parameter C/{(i). 

Therefore the total contribution of Qi x Qi to this term does not exceed 
Cmax(zi, zIbIl^). 

On the other hand, contribution to the term in question of each couple 
of balls contained in Q.2 x Q2 does not exceed its contribution to (5.3.12) 
and therefore the total contribution of Q2 x Q2 to this term does not exceed 
expression (5.4.13). 

So, term (5.3.11) with A = A/v does not exceed (5.4.13). 

Therefore we arrive immediately to 

Theorem 5.4.4. Let M > 2, B < and condition (0.4.1) be fulfilled. 

Then under assumption (5.4.1) 

(i) The following estimate holds: 

(5.4.14) E^^ < S^"" + (Tr{HA,w - u)' + J Pell/l/""" + i^) dx^ + 

Cmax(zi, ziei/.^); 

(a) As a > Z~^ The following estimate holds: 

(5.4.15) E"^"" <£^^ + Scott + C max (Zi, zifii/.^) + CZ^fii + Ca^^zi; 

as a < Z"^ one should replace the last term in the right-hand expression by 
CZ^ and skip Scott; 
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(in) As B < Z and a > Z I 

(5.4.16) E^^ < 8'^^ + Scott + Dirac + Schwinger+ 

CZi (Z-^ + [BZ-^Y + {aZ\)-^) . 

Here proof of (iii) is due to the same arguments as in the case 6 = 0. 
Combining with the estimate from below we also conclude that 

Theorem 5.4.5. (i) In the framework of theorem 5.4-4 the following esti- 
mate holds: 

(5.4.17) D(py;-p"^^ Pv-P^^) < Cmax(zi, zifii/.*); 

(a) In the framework of theorem 5.4 -4 (Hi) the following estimate holds: 

(5.4.18) Dip^-p"', p^-p''') < CZl{Z-' + {BZ-y + {aZ'3)-'). 

5.4.3 Estimate for D-terms under assumption (5.4.3) 

Let assumption (5.4.3) be fulfilled. Let W = W^, and i = £,yhe a potential 
and a scaling function (used to derive semiclassical remainder estimates) 
for this u < (and N < Z) while Wq and Iq be a potential and a scaling 
function for z/ = (and N = Z). 

Let us start from the trivial arguments. Note that 

(5.4.19) |P^(1/1/ + A) - P'b{W + v)\ < CW'^\-v\ei + CB\\-v\^e2 

where Ox and 82 are characteristic functions of = {x : W{x) >\v\} 
and {x : < l/l/(x) < respectively. Let 

A D ([P^( 1/1/ + A„) - 1/1/ + u)\e,, [P^( 1/1/ + A/v) - 1/1/ + vW,) . 

Then in virtue of (5.4.19) 

Ji < CD(l/l/^^i, l/l/^^i)|A/v - 

Note that D(l/l/5^i, l/l/^^i) x ((Z - A/);^ 6"^, Zifi-f) in our cases and 
using (5.4.4) one can prove easily that Ji < CZs. 
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However estimate for a contribution of zone {x : < 1/1/ (x) < Co|i^|} is 
much worse: 

(5.4.20) J2 < CB^D{e2, ^2)|A/v - i^l < CB^r^{\u\/ B^)hl\XN - 

1 

where for B < [Z — N)\_ wc should replace (|z/|/6^)2/.^ by r^. Using (5.4.4) 
we conclude that for (Z - N)l < B < 

J2 < CerV|t/.2max(zi, e^) x Ce(Z - /V)f m max(zi,S5/.i)/.^, 

and therefore we arrive to the last two cases below; the first case is proven 
similarly: 

({z-n)-JzIb\ 

(5.4.21) J2<cl (z - /V)f max(zi, Bhi)Ll 

^ [z - n)Iz^obIlI 

This is not a good estimate. To improve it let us note that 

(5.4.22) If |Aa/ -iy\< C max(ei, (Z - A/)f ) then J2 does not exceed (4.4.6) 
and therefore we can assume that 

(5.4.23) lAw-i^l > Cmax(6i,(Z-/V)f). 

Proposition 5.4.6. (i) Let (Z - N)l < B < Z^ and 

(5.4.24) CoSi < |A/v - < QB^^l'^, 
then the truncated semiclassical error does not exceed 

(5.4.25) M := CZi + CB?^{\i^\/ B^)'^ L x {B-^\Xn - 

(a) Let iZ-N)l<B < Z^ and 

(5.4.26) Cifi^|i/|? < |A/v-i/| < CiS5|i/|U, 
then the truncated semiclassical error does not exceed 

(5.4.27) M CZi + CS?^/.; 
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(in) Let B < {Z - N)l and 



(5.4.28) 



(Z - N)l < \X„ -u\< Ci(Z - N) 



then the truncated semiclassical error does not exceed 



(5.4.29) 



M := CZi + C{Z-N)1\\n-v 



-1. 



4 



(iv) Let B <{Z- N)l and 



(5.4.30) 



Co(Z - N)+ < |A/v - i/| < CiZi(Z - N)l 



2 



then the truncated semiclassical error does not exceed M := CZs. 
Proof. Easy proof using arguments of the proof of proposition 5.4.1 is left 



Proposition 5.4.7. In the framework of proposition (i)-(iv) term (5.3.12) 
does not exceed CM^{B\\f^ — v\2)3 + (4.4.6) with M defined defined in the 
corresponding cases in Proposition 5.4-6. 

Proof. Using proposition 5.4.1 one can prove easily that 

(5.4.31) The contribution of {x : mirim \x — y^l < er} to (5.3.12) does not 

5 

exceed CZs. 

Now we need to estimate excess of expression (5.3.12) over semiclassical 
D-term (with X = u) which has been estimated by (4.4.6). To do so we need 
to estimate 



(5.4.32) D{[Pb{W + u) -Pb{W + \)]e, [Pb{W + u) - Pb{W + A)]^) 



which is the contribution of the domain Q' := {x : £{x) < CB -^{X — u\} 
where 6 is the characteristic function of Q'. Recall that on the complimentary 



domain |Pe(l/l/ + z/) - Pe(l/1/ + A)| < CE''^. Let us consider 

(5.4.33) D([Pe(l/l/ + z/)-Pe(H/ + A)]^o, [Pb{W + u) - Pb{W + X)]eo) 



to the reader. 



□ 



85 



and 

(5.4.34) D{[Pb{W + u)- PB{W + X)]9t, [Pb{W + u) - Pb{W + X)]et') 

where is a characteristic function of 

Q() :={x: e{x) < to := {\Xn - iy\B-^)^ 

and 6t is a characteristic function of Q.[ := {x : t < i{x) < 2t} with 
t>t'> to. 

When calculating (5.2.5) the contribution of Qg is x 6|A/v — i^|2 mes(Qo) 
and therefore due to proposition 5.4.6 

(5.4.35) mes(Qo) < CM{B\Xn - 
while term (5.4.33) is 

X fi2|A/v-i/|D(^o,^o) < Cfi^|Aw-i^|(mes(f2o))^ < C/wi (S|A/v - z/j^) ^ 
where the middle inequality 

(5.4.36) D(xg,Xg)< C(mes(G))^ 

is well known^^) and the last one is due to (5.4.35); xg denotes characteristic 
function of G. 

Similarly, when calculating expression (5.2.5) one can see easily that the 
contribution of Q'j. is x 6|A/v — u\{B'^t^)''^ mes(Q() and therefore 

(5.4.37) mes(f2;) < CM|Aa, - 

while term (5.4.34) is x |A/v — i^|^t'~^t'~^D(^t, 9t') which does not exceed 

(5.4.38) C|A/v - uft-^t'-^ mes(Qt) mes(fitO [ 
due to inequality 

(5.4.39) D(xG, Xc) < C mes(G) mes(G') [max(mes(G), mes(G'))] 

Really, among uniform solids of equal mass and density the ball has the least 
potential energy; then C = |(127r)3. 
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which trivially follows from the obvious inequality D{xg, 5z) < C(mes(G))3, 
where 5^(x) = 6(x — z). 

5 2 

Due to (5.4.37) expression (5.4.38) does not exceed CMsf^s; recall that 
t > t'. Since summation with respect to f > t' and then with respect to 

5 _2 

t' > to returns CMit^ ^ we conclude that term (5.4.33) with replaced 
by 9" (the characteristic function of {x : £(x) > to}) also does not exceed 
CMI{B\Xn-u\'2)i. □ 

So, we have now two estimates for an excess of expression (5.3.12) over 

5 11 

(4.4.6): one estimate is CMs (S|A/v - u\2)3 with M = /W(|Aa/ - derived 
in proposition 5.4.6 and another one is due to (5.4.20). Let us consider the 
best of them. Note that the former estimate consists of two terms each due 
to the corresponding term in the definition of M. The second term in the 
framework of proposition 5.4.6(i) is 

C(fi5?^|i/|?/.|A/v - [B\Xi^j — X B^?^\u\^L^X|\J — 

and taking minimum of this expression and CBr^lul^ L'^\Xn — i/| we see that 
this minimum does not exceed 

(5.4.40) C(eTrT|zy|i5/.i)i(erV|5/.2)^ x CBh' {Z - N)^ L^s x 

^ f (z-/v)|e/.i (z- /v)| < s < z|, 

\{Z - N)^zf'Bh?^ Zl<B<Z^ 

11 _ 2 I I 1 2 

which is achieved as \Xn — x Sisris |i/|~3o/.~i5 and one can see easily that 
it is in the framework of assumption of proposition 5.4.6(i). 

2 5 11 

Furthermore, the first term is {Z3)i{B\XN — z/|2)3 and it reaches its 

2 1 1 

maximum at the largest value of |A/v — (which is max(Z3, B2L)\v\^) and 
we leave to the reader to check that it does not exceed (5.4.40). 

In the framework of proposition 5.4.6(iii) one should replace (z//6^)^ by 
? and L by 1, so B2j'^\v\^ L\Xn — B^?^\Xn — further, one should 

preserve S|Aa/ — z^l^ and therefore the second term becomes 

{B^r^\XN - {B\Xn - x B'^-r^\XN - 

which should be run versus (6|A/v — -\- \v\^\Xn — v\f'?^ resulting in 

CBh^ + Cfi5i|i/|7 7=5 ^ cB^Z - /V)7 + CS5i(Z - /v);^ 
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and both of them are less than in (4.4.6). Meanwhile the first term becomes 
ZtS^|A/v — < {Z — N)^Zt B^. Combining with (5.4.40) we arrive to 



(5.4.41) 



A _ii _i 



S < (Z - N)l, 
{Z-N)i<B<zl 



Zl<B<Z^ 



{Z-N)^ZtB-3 
{Z-N)^BLr5 
^ {Z - N)^ zf' Bh'i 

6AA Summary 

Theorem 5.4.8. Let M > 2, B < Z^. Then 
(i) The following estimate holds: 

(5.4.42) E^^ < S^"" + (jr{HA,w - 1^)' + J Pe(l/1/^'' + i^) dx) + CR+ 

(4.4.6) + (5.4.41); 

(a) The following estimate holds: 

(5.4.43) E^^ < + Scott + C/?+ CZt si + (4.4.6) + (5.4.41); 

(Hi) As B < Z and a > Z~i 

(5.4.44) E^^ > S'^^ + Scott + Dirac + Schwinger+ 

CZi(Z-'5 + (eZ-^)^ + (aZ5)-^) 

Remark 5.4-9. Unfortunately there are gaps in the proofs of V. Ivrii [Ivr2] 
in the case of M > 2 and large Z — N > and I was unable to fill them. 
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6 Negatively charged systems (coming) 



7 Positively charged systems (coming) 



8 Appendices 

8. A Electrostatic inequalities 

There are two kinds of electrostatic inequalities: those which hold for 
any fermionic state U' and those which hold only for ground-state (or near 
ground state) Inequalities of the first kind do not depend on the quantum 
Hamiltonian and they are (24.2.1) of [Ivrll] repeated here: 



(8.A.1) Y 



l<j<k<N 





and (8. A. 5) below. 

Inequahties of the second kind e = 0: 



(8.A.2) Yl 



l<j<k<N 




-D(pv|/,pvi/) - CZ3, 



and more precise one (8. A. 26) below 

As 6 = const there is inequality established in E. Lieb, J. P. 
J. Yngvarsson [LSY2] (p. 122): 



Solovej and 



Theorem 8.A.I. Let B = const. Then for the ground state U' 



(8.A.3) 



/ 



pldx< CZiA/^(Z + A/)3(l + 
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In particular for c < Z < cN the right-hand expression does not 
exceed 



(8.A.4) CZt(l + eZ-^)i X C 



Zi as e < Z5, 

Zisgs as B> Z5. 



We want to establish inequality similar to (24. A. 2) of [Ivrll] but in 
magnetic case. We will use for this the following 

Theorem 8.A.2. ^4) Fix<d<6 < 1/6. Then for any density matrix G an 
any density po{x) > we have 

(8.A.5) Yl [ \XJ - Xkl^'l'^ixi XA/)|'C/Xi---C/XA/ > 

l<j<k<N •' 

D(po,P7) - ^D(po,Po) - // \G{x,^;y,^')\^\x-y\'^ dxdy 



-C\\pfj;-\\p\\T^'v{^,Gfl^'^ 
where p = po + pc + Pw, v{^, G) := Tr{^{l - G)), 

(8.A.6) 7 = 7y,(x,y) = A/ J ^{x,X2, ... ,XN)^\y,X2, ... ,xn) dx2- ■ ■ x^ 

is two-point one particle density. 

Recall that ||.||p denotes ^''-norm. 

Connection between (8.A.1) and (8. A. 5): as we set G = we get ^ = 
and the last term in (8. A. 5) becomes IIPII5/3 • IIpHi^^; on the other hand, 

IIPII4/3 < IIpIIs/s ■ llpliy^) so (8. A. 5) is slightly deteriorated (8.A.1) with G = 
but with "free" po. 

Let us follow G. Graf and J. P. Solovej [GS] further albeit in the case of 
magnetic field. Let us estimate first ||p||5/3. 

24) Lemma 6 of G. Graf and J. P. Solovej [GS]. 
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If p = p^^ direct calculations show that as N >^ Z 

(8.A.7) J p^^ dx^m\n{Z,N), 

(8.A.8) j (p"''')5c/xx Cp*tr*3 = czi(l + eZ-5)^ 

(8.A.9) y" £/x X Cp*ir*3 = CZl{l + BZ-'^y 
with 

(8.A.10) r* = min(Z-^ e-izi) x Z"^ (l + eZ-s)"^ 

(8.A.11) p* = min(A/,Z)r*'3 

and we have \\p\\l^/l ■ \\p\\y^ x \\p\ff^l for p = 

If p = Pm/ we use magnetic Lieb-Thirring inequality (see f. e. L. Erds [E]) 

(8.A.12) MHa.w)' J PB{W)dx 

and therefore 

(8.A.13) {m,^\f) >Tr- {Ha, w) + J WpM,dx 

- J Vp^ dx + ^D{p^^,, pm,) - C||pxu||J^3 
which due to (8.A.12) is greater than 

(8. A. 14) J {-CPb{W) + l/l/py,) dx- 

J VpM, dx + ^D(py,, pm,) - CllpvllJ/3 > 

3eo J TB{p^)dx- j Vpx„c/x+ ^D(pvj;,pxi,) - C||pvj/ 114^3. 

as we picked up 1/1/ : CP'g{W) = px^. 

The first two terms in the right-hand expression are estimated from 
below by 

2eo J TB{pM,)dx-C J PB{V)<j)dx-C J VpM,{l-<P)dx 
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where supp(0) C {x,i{x) < 2r*} and supp(l - (b) C {x,i{x) > r*}. 

One can see easily that the absolute value of the second term is x 

7 / 4 \ - 

Za (^1 + BZ^ij ^ while the absolute value of the third term does not exceed 
CZ f V(l — 0) dx X CZ^r*~^ which does not exceed the same expression 

Zi{l + fiZ-t)^. Therefore 

(8.A.15) (Hv|/, \|/) + CiZi (1 + BZ-ly > 

2eo J tb{pxv) dx + ^D(px„, p^,) - C\\pM,\ff^l. 

Note that ||pvu|l4/3 calculated over domain {x : pv];(x) > B^} does not exceed 
^11/^^ II 5/3 ■ II pill with norms calculated over the same domain which does 
not exceed CT^Z^ with T = J tb{p^) dx. 

Meanwhile ||pv]y||4^3 calculated over domain {x : p\v{x) > B^} does not 

1 5 

exceed C||p||| ■ with norms calculated over the same domain which 
does not exceed CZ^fis Te. 
Therefore 

(8.A.16) l|pvi;||J^3 < cr^z^ + czie^ 

and therefore (8.A.15) imphes that if 

(8.A.17) (HM/, M/) < QZl (l + BZ''^) ^ 

then 

(8. A. 18) T = J TB{p^v) dx < C2ZI (1 + eZ-t) i 
and 

(8.A.19) \\pXf,<CZ'/\l + BZ-ly; 

taking G = we arrive to (8. A. 3) as A/ x Z. 

However on our preparatory step we need to estimate also UpvuHs/s and 
due to (8. A. 18) we need to consider only norms over {x : p\\;{x) > Ssj. Then 
II 1 1 5/3 < ^I|py|l4/3"'^^||pvi/||3^^ < and plugging the same estimates (8. A. 19), 
(8. A. 19) we conclude that 

(8.A.20) ||pM,||5/3 < CZ^/3(l + SZ-5)^ 
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Now we assume that B < Z^, take G = e(x,y, z/) where e(x, y, i/) is 
the Schwartz kernel of spectral projector for potential W approximating 
1/1/"'"'^ and z/ < is a chemical potential. One can prove easily that IIPGII5/3 

satisfies the same estimate and we need to estimate Tr(7y/(/ — E(/x))). 
Consider 

(8.A.21) N{Ha,w(.,)^, ^) - MHE{iy)) - aTr(7^(/ - E{u))) 

>f Pd^TrEiP)- f {P-u + a)d0TrEiP) 

J/3<0 J P<v 

= - I (/9 - ^ + a)dpE{^) 

J u—a<l3<u 

= -aE{i,) + f E(/3) d^. 

J v—a<fi<v 

We can replace E{j3) by / P'{W + j3) dx with a resulting error 0{Zafi^), 
h := BZ~^. Then the right-hand expression becomes 

(8.A.22) -L{a):= J (^-aP'^iW + u) + P'b{W + u - /3) d dx = 

-j\a-P)(^j P'^{W + u - ^) dx) 

Therefore 

(8.A.23) a(Tr{-f^{l - E(/x))) - CZH^) < 

N{Ha,w(.,)^, ^) - MHE{iy)) +L{a). 

Note that adding to the selected terms — |D(p"'"'^, p^^) we obtain exactly the 
snippet occurring in the lower estimate of E/v but in virtue of the upper 
estimate it should not exceed Q = CZs (l + 6Z~3)5 < CZs + Ch^Zi and 

therefore plugging a = Z^H^ we conclude that 
(8.A.24) Tr(7xu(/ - £(/x)) < ZH^ 

provided we prove that 

(8.A.25) L{a) <Q as a = zh^. 

Therefore modulo proof of (8. A. 25) we arrive to the estimate (8. A. 26) below: 
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Theorem 8. A. 3. Let N ^ Z and B < Z. Then for the ground state energy 
(8.A.26) Yl /kj-^/cr'l^(^i.-.^/v)l'c^^i---c/^A/ > 

l<j<k<N 

^D(pT^ p""") + Dirac - CZl-\l + B^) 

To prove (8. A. 25) we note that 0<Pg(i/i/)xi/i/5-|- 6m/^i. One can prove 

7 13 

easily then that L{a) < Ca'^ + C64q;2 which obviously implies (8. A. 25). 

8.B Very strong magnetic field case 

Let us consider now the case Z^ < B < Z^. 

Proposition 8.B.I. Consider Schrddinger operator Ha,w with a constant 
magnetic field of intensity B and potential W: W < Z\x\~^. Let 0(x) := 
(f)r{x) be r -admissible function. Then as Z^ < B < Z^ and r x Z~^ 

(8.B.1) \e{x,y,0)\ < CZB inB{0,r) 

and 

(8.B.2) All eigenvalues are > —CZ^. 

Proof. Without any loss of the generality one can assume that 

(8.B.3) Ha,w = Dl + Dl + (Di - Bx2f - W. 

Consider f e sL?] then ||E(A)f || < \\f\\ and then one can prove easily (8.B.2) 
and that 

(8.B.4) \\HA,oE{\)f\\ < {CZ^ + \+)\\f\\. 

Really, + CZ^ > 1/1/ in the operator sense. 

Then (8.B.4) implies that in 6(0, r) x B{y' , r') with r' = 

|P"E(A)f I < CZ"3S^I"'I Va : |a| < 2 VA < 

with P = {Dl - Bx2, D2, D3) and therefore ||E(A)/'||^ < CZifiiH/^H. Then 

\\E{x, .,X)\\^2\\ < CZsSs. 

Repeating the same arguments with respect to y we arrive to (8.B.l).n 
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Corolleiry 8.B.2. In the framework of Proposition 8.B.1 with e =5f°°(6(0, r)), 

\^oo < 1 



(8.B.5) I J <l){x)e{x,x,0)dx\ < CZ^^B, 

(8.B.6) D(0(x)e(x, x, 0), 0(x)e(x, x, 0)) < CZ'^fi^ 

and 



(8.B.7) \f f (f){x)e{x,x,T)dTdx\< 

J —oo J 



CB. 



8.C Riemann sums and integrals 

life ^°°(R+) and fast decays at +oo then 

POO 

(8.C.1) /^(0)/j + ^2/^(2n/j)/j- / f{t)dt + Y,f^mf^^'"~^\^)h^'"' 

n>l m>l 

f{t) cyt + ^<f(2'"-^)(o)/j2'" 

n^u - m>l 

as /) — )■ +0. The proofs of both formulae foUow from Taylor decomposition 
and observation that the odd powers of h should disappear. Taking f{t) — 
e"^^/^ with Rez > we arrive to 

(8.C.3) l-^-~E 

sinhlz) ^-^ 



m>l 

2m 

sinh(z) 



m>l 

as |z| <S 1. In particular, = | and k'^ = — | 
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